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ABSTRACT 


This research examines the potential for the conjunction between photovoltaic 
(PV) and thermophotovoltaic (TPV) technologies for spacecraft power production. There 
is sufficient overlap between the sources of energy used for these devices and the 
function of the devices themselves that either PVs or TPVs could gain improvements in 
efficiency from the integration of the other type of device, or that a hybrid device could 
be developed. As a proof of concept, a GaAs PV cell and GaSb TPV cell were modeled 
in a tandem design using Silvaco ATLAS, with varying PV cell substrate thicknesses, 
and simulated under the AMO spectrum to determine the potential range of efficiency 
gains for a PV device integrated with a TPV device. The same design was then 
tested under a 2000 K blackbody spectrum—to approximate use in a radioisotope 
thennoelectric generator (RTG)—to detennine if similar efficiency gains could be seen 
for a TPV device integrated with a PV device. The possible gains with a PV-TPV design 
under AMO are clear, potentially resulting in cells with a 30-34% overall efficiency. The 
possible gains for a PV-TPV device utilizing a blackbody spectrum are less clear, and 
would benefit from further design and investigation. 
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I. 


INTRODUCTION 


A. BACKGROUND 

Modern spacecraft power production is dominated by the use of photovoltaic 
(PV) arrays, which draw electrical energy from the most prominent power source in our 
solar system, the Sun. These arrays are large, and pose considerable issues with 
spacecraft control and maneuvering; the larger they are, the more power is required for 
the spacecraft to maneuver and maintain attitude control. Any potential efficiency gain 
for these arrays, even as little as 1%, results in significant decrease in the overall array 
size required, which leads to lower power requirements for maneuvering and attitude 
control, which in turn leads to the potential for further reduced array size. As the 
maximum mass that can be sent to any given orbit is limited by current space launch 
technologies, lower mass requirements in any given spacecraft system, particularly the 
power system, mean the potential for increased payload mass and functionality. In turn, 
higher efficiencies in the same size of array also mean more power available for the 
spacecraft payload, while requiring no or little additional mass, translating to similar 
potential increases in payload functionality. Thus it can be seen that even a small 
electrical efficiency increase in PV cell design can result in a cascading improvement to 
the overall quality or functional potential of a spacecraft. 

It is for this reason that various PV cell manufacturers continue to seek even the 
most minor improvements in efficiency in even the most advanced PV cell designs. 
Various PV cell technologies have improved in efficiency by as much as 1% from year to 
year, but many have taken decades to improve a similar amount, and some seem to have 
leveled off entirely [1]. Increasingly, PV cell designers are looking for ways to merge 
various technologies or go outside the normal PV design scheme to enhance power 
production. One potential improvement, associated with developments such as tandem, 
dual and multi-junction designs, is the use of semiconductor devices with band gaps of a 
magnitude such that longer wavelengths than PV cells are normally designed to utilize 
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can be absorbed. Devices that utilize this area of the electromagnetic (EM) spectrum for 
power production are kn own as thenno-photovoltaics (TPV) and are not only a relatively 
new area of exploration in power production, but the use of some materials for TP Vs is 
largely unexplored relative to their use in conjunction with PV materials. 

One area where the use of TP Vs is being explored vigorously is in relation to 
radioisotope thermoelectric generators (RTGs), which use heat generated by the decay of 
radioactive materials to produce electrical power, and are the second most common 
device used for spacecraft power production. The devices previously utilized for power 
conversion in these systems, thermo-electric generators (TEGs), require a fairly extreme 
temperature gradient to operate effectively. Alternatively, the blackbody radiation 
emitted in these devices can be converted to electricity via TP Vs, which have a much 
higher potential efficiency. 

The nature of the majority of photonic energy sources from which either PVs or 
TP Vs draw power is such that the distribution of energy emitted by the source possesses 
a dominant wavelength, but shows energy distribution across the range of the EM 
spectrum similar to a nonnal distribution skewed toward that frequency. This is 
commonly referred to as a blackbody curve, and this distribution of energy across a 
black-body curve means that there is a decent quantity of energy distributed across the 
EM spectrum for any given energy source. A more flexibly designed device employing 
both PV and TPV technology would not only serve to generate a greater quantity of 
energy from a given source, but could also be more flexible in the types of sources from 
which it could generate electricity. Given the consistent drive for improved efficiency in 
PV power production, and the need for flexible or adaptable power production systems, it 
would seem clear that this possibility is worthy of investigation. 

B. PROBLEM STATEMENT 

The commonality between how TPV and PV arrays operate, as well as the 
potential for improvement for each type of device, involves overall broader access to 
areas of the EM spectrum normally considered separately or specific to one type of 
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device or the other, such as visible light being used for PVs and infrared (IR) radiation 
emitted by the decay of a radioisotope in an RTG being used for TPVs. The fact is that 
both these areas of the EM spectrum and these technologies have considerable overlap in 
research and development, and considerable potential for enhancing power production in 
either type of system. The key to this may simply be the integration of PV and TPV 
devices into a single PV-TPV cell. 

The means for accomplishing this are fairly straightforward, given sufficient 
knowledge of semiconductor device manufacturing. It essentially involves the layering of 
one semiconductor device material on top of the other. While there are a variety of 
manufacturing techniques that could be involved in accomplishing this, they all typically 
result in one of two major types of design. The first of these is a stacked or tandem 
design, in which a device composed of one material type is covered by a transparent 
insulating layer, for electrical isolation, and the other device is then laid on top. This 
allows power to be generated by the upper device, and then for any EM radiation to 
which it is transparent to pass through to the lower device. The second major design type 
is a dual or multi-junction design. This design is essentially a more modern extension of 
the tandem design in which the semiconductor material layers of one type are deposited 
or fonned directly on top of the other type; though not electrically isolated, per se, the 
layers are typically separated by two thin intermediary layers of opposite polarity which 
form a tunnel junction between the two, allowing electrical current to pass through from 
the upper layers to the lower layers in only one direction. 

Given the proper manufacturing facilities, multi-junction designs are typically less 
expensive to produce than tandem designs; however, their modeling is typically more 
complex and does not allow for the independent assessment of individual layer or device 
performance, which can make optimization more difficult. Also, while it may be more 
expensive to manufacture a tandem design, modeling a tandem design vice a dual¬ 
junction design poses no additional cost, and the potential output from a design of a 
tandem PV, TPV, or PV-TPV cell can be mathematically associated with the potential 
output for a multi-junction design—where the same cannot be said of a multi-junction 
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model. To this end, while it might be the more expensive option to manufacture, it would 
seem that use of a tandem design model would result in the greatest flexibility for overall 
technology assessment. 

C. METHODOLOGY 

For the assessment of the potential gains in efficiency and power production to be 
had by the conjunction of PV and TPV power technologies, a tandem PV-TPV device, in 
which a PV cell is positioned above a TPV cell, was modeled parametrically using the 
Silvaco ATLAS technology computer aided design (TCAD) suite to determine the range 
of potential power output. The material chosen for the PV cell was gallium arsenide 
(GaAs), which is fairly commonly employed in the PV industry to produce PV cells of 
relatively high efficiency. The material chosen for the TPV cell was gallium antimonide 
(GaSb), which has fairly well known and established semiconductor material properties, 
but has not been widely employed or investigated as a TPV material or as a potential PV 
material. The assessment of devices modeled using these materials should provide a solid 
assessment of the potential of PV-TPV devices, as well as a stable overall base from 
which to draw comparison of future devices employing other materials. 

PV and TPV devices were modeled and assessed independently under both the 
AMI.5 and AMO spectrums, for validation of the expected output from these models. The 
devices were then combined into a tandem model for assessment of overall device 
potential, and various device characteristics such as doping and substrate thickness were 
varied to determine a range of potential device efficiency and power production, as well 
as provide a rough approximation of optimal device parameters and an overall assessment 
of the potential gains from integration of PV and TPV devices. The results of that 
research are the subject of this thesis, and are provided as follows. 

Chapters II and III of this thesis discuss the nature of solar and radioactive energy 
sources and the means by which the photonic energy from those sources is directly 
converted to electrical energy, specifically via PVs and TP Vs, with some commentary on 
TEGs—as they are employed in RTGs—and how their performance compares to TPVs. 
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Chapter IV of this thesis details the modeling software, designs, and assumptions 
employed in evaluating the potential of PV, TPV, and tandem PV-TPV devices. Chapter 
V discusses the results of these evaluations, and Chapter VI discusses conclusions drawn 
from these evaluations and provides recommendations for future work in this area. 
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II. ENERGY SOURCES 


A. SOLAR RADIATION 

The predominant source of energy in our solar system is the sun. In conjunction 
with the Earth’s rotation and the resulting solar gradient, energy from the sun drives 
thennal atmospheric and ocean currents around the globe. Wind and ocean wave 
mechanisms have been developed to access the energy that flows through these currents. 
The cycle of evaporation, atmospheric condensation, precipitation, and water flow across 
land—commonly referred to as the water cycle—is also driven by solar radiation. 
Hydroelectric power has been developed to convert the potential energy of water, as it 
accumulates and cascades through a drainage basin, into kinetic and electric energy. Even 
fossil fuels are a tertiary conversion of solar energy, having been converted by plants 
through photosynthesis to carbohydrates and cellulose, sometimes by animals into more 
carbohydrates and protein, and then being geologically converted by heat, pressure, and 
time into vast pools and bubbles of hydrocarbons. A bevy of combustion mechanisms and 
engines, of course, allows us to access this stored energy. One of the most direct methods 
for converting solar radiation to usable energy, however, is photovoltaics (PV). 

The photoelectric effect, the excitement of electrons in a material by incident 
light, was first discovered in 1839 by Edmund Becquerel [2]. Over the past one-hundred- 
and-seventy-six years, various endeavors in research and development have transformed 
this observation into an industry for the design and manufacture of crystalline materials 
that directly convert the electromagnetic radiation of the sun into electricity. 

Terrestrial use of PVs for power production is becoming remarkably widespread, 
but is still limited in effectiveness by such variables as atmospheric losses of solar 
irradiance, the terrain, weather, and latitude of the installation, as well as the seasonal 
cycle and time of day. Figure 1 describes the overall annual impact of some of these 
variables on direct solar irradiation. The resulting lack of consistency in solar irradiance 
at any given terrestrial location results in the construction of larger arrays, more extensive 
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solar energy production facilities, and considerable energy storage facilities to regulate 
the power available throughout the day, as well as limiting the sites that appear viable for 
employment of solar technology. 


Figure 1. World Map of Direct Normal Solar Irradiation (Insolation) 
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This map describes the considerable variation in annual solar insolation throughout the 
terrestrial environment. Provided by SolarGIS (solar geographic infonnation system) 
http://solargis.info/doc/_pics/freemaps/1000px/dni/SolarGIS-Solar-map-DNI-World- 
map-en.png 

By contrast, spacecraft PV power production systems are not so limited. With the 
exception of spacecraft in low-Earth orbit (LEO)—which frequently pass through the 
Earth’s shadow and have an approximately 60 minute-30 minute day-night cycle—the 
majority of spacecraft using PV systems, specifically those in geosynchronous, 
geostationary, and medium earth orbits (GSO, GEO, and MEO, respectively), are able 
to maintain almost constant power production. Also, the lack of atmospheric losses 
means that the solar irradiance received by a spacecraft PV system is approximately 
30% greater, overall, than that of a terrestrial PV system. The only real limiting factor for 
a spacecraft PV system is the relationship among the needs of the spacecraft, the lift 
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capabilities of modem launch systems, and the efficiency of the PV cells used— 
discussed in depth later—which is the main driver in the overall efficiency of the PV 
system. 

B. RADIOACTIVE DECAY 

The second most abundant source of energy in our solar system is heat generated 
by the decay of radioactive materials. Thirty-eight of the one-hundred-and-eighteen 
elements on the periodic table are radioactive, and a recent study by the KamLAND 
Collaboration [3], which employed geo-neutrino flux to determine a more accurate model 
of the Earth’s core, revealed that decay of uranium-238 and thorium-232 may contribute 
to almost half of the Earth’s geothennal heat flux (approximately 20.0 of 44.2 TW). 
Though this obviously contributes to our use of geothennal energy for heat and electricity 
production, more direct methods of converting radioactive heat to usable energy have 
been employed for decades. For example, nuclear power plants use a controlled nuclear 
reaction, an accelerated form of nuclear decay, to heat water to steam to run a turbine 
driving an electrical generator. However, direct conversion of heat from spontaneous 
radioactive decay to electrical energy is the province of thermo-electric generators and 
thenno-photovoltaics. Both of these technologies have been examined in conjunction 
with an array of radioactive elements for use in energy production [4], [5], [6], [7], 

The process of radioactive decay has simply to do with the breakdown of unstable 
atoms or unstable isotopes of an atom, through the emission of one of three types of 
radiation, and the formation of more stable atoms or isotopes thereafter. The three types 
of radioactive emissions are alpha particles, which are essentially a helium-4 nucleus, 
beta particles, which are emitted electrons or positrons, and gamma rays, which are 
extremely high energy photons. Each of these is of progressively higher concern from a 
safety perspective, as their lower mass or higher energy results in greater potential 
damage to surrounding materials. As a result, shielding, involving thick layers of dense 
materials such as lead, through which the radiation will be less likely to penetrate, is 
typically employed to absorb emitted radiation and protect surrounding personnel and 
equipment. Electronics are particularly sensitive to radioactive emissions, as the emitted 
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high energy particles and electrons can cause either permanent damage to micro or 
nanoscale electronic components or simply cause an upset, or error, in the status or 
processing of an electronic device. In fact, the abundance of such radiation in the space 
environment is one of the justifications for the high cost of space-qualified components. 
These components have to be sufficiently radiation shielded or redundant to survive for 
the life expectancy of the spacecraft, which can be as much as ten to twenty years. 

Despite the apparent concerns over utilizing radioactive materials as a thermal 
heat source, the practice is not uncommon, particularly in environments or applications 
where the risk to human health or life is minimal. Some of the more prominent of these 
applications are as the power source for NASA deep-space craft [8], These craft are 
launched for missions sufficiently far from the sun that the intensity of its radiation is 
significantly diminished, which means a photovoltaic array would need to be ridiculously 
large relative to the size of the spacecraft in order to provide sufficient power. Likewise, 
consumable fossil fuel or hydrocarbon sources are of sufficiently low energy density that 
they cost almost as much, if not more, energy to put into space than they would provide 
for the space craft, and the volume that could be launched with the craft would not be 
sufficient to supply it for the duration of its mission. For example, NASA’s recently 
popularized New Horizons mission [9], to Pluto and the Kuiper Belt, was launched in 
January of 2006, nearly ten years ago, and is expected to operate for at least another five 
years (through 2020), which is not uncommon for a NASA deep-space mission timeline. 
The energy density, long life, and consistent energy production of a radioactive heat 
source make such missions possible. 

Radioactive thermal heat sources have even been used by the United States Navy 
(USN), and investigated for use in forward operating applications by the United States 
Marine Corps (USMC). Table 1 compares general physical characteristics for various 
radioactive materials, and Table 2 shows mass and cost comparison statistics for some of 
the same radioactive elements, the analysis of which was conducted by Langham [4] at 
NPS in 2013 for a prospective 300 W RTG power supply. It can be seen that some of 
these elements are more highly available, less expensive, provide a greater power output, 
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or have a half-life more or less conducive to space applications. Pu-238, one of the 
radioactive isotopes of Plutonium (Pu), has some of the better characteristics for a long 
tenn power supply, and is, in fact the fuel normally used in a NASA deep-space RTG. 


Table 1. Radioisotope Data 


Isotope 

Particle 

Emitted 

Specific 
Power (W/g) 

Density 

(g/cm3) 

Compound Specific 
Power (W/g) 

Half-life 

(yr) 

Cost 

($/g) 

Availability 

(kg/lOOMweY) 

Sr-90 

Beta 

0.96 

2.6 

0.22 

27.7 

20 

16 

Cs-137 

Beta/Gamma 

0.42 

3.42 

0.12 

30 

6.5 

36 

Po-210 

Alpha 

141.3 

9.298 

134 

0.379 

2800 

1 

Pu-238 

Alpha 

0.56 

19.86 

0.39 

87.48 

300 

15 

Cm-242 

Alpha 

120 

13.51 

98 

0.45 

2000 

1 

Cm-244 

Alpha 

2.65 

13.51 

2.27 

18.1 

170 

1 

Tm-170 

2e:a 

13.6 

9.32 

1.2 

0.35 

136 

1 

Ru-106 

Beta 

33.1 

12.45 

1.1 

1 

120 

2.3 

Pm-147 

Beta 

0.33 

7.26 

0.27 

2.6 

75 

3.6 


This description of characteristics of general concern in relation to a variety of 
radioactive isotopes is from R. C. Langham. “Feasibility study and system architecture of 
radioisotope thermoelectric generation power systems for USMC forward operating 
bases.,” M.S. thesis, Dept. Sys. Eng., Naval Postgraduate School, Monterey, CA, 2013 


Table 2. Initial Cost and Weight (for a 300 W RTG) 


Isotope 

Pure isotope 
mass per 
300W (g) 

Compound 
mass per 
300W (g) 

Estimated 
cost per 
300W (S) 

Initial pure 
mass for 300W 
at EOL (g) 

Initial compound 
mass for 300W at 
EOL (g) 

Initial estimated 

cost for 300W at 
EOL (S) 

Sr-90 

313 

1364 

6250 

401 

1751 

8027 

Cs-137 

714 

2500 

4642 

900 

3150 

5850 

Po-210 

2 

2 

5945 

186090003 

196227741 

521052007084 

Pu-238 

536 

769 

160714 

580 

833 

173967 

Cm-242 

3 

3 

5000 

12231943 

14977890 

24463887084 

Cm-244 

113 

132 

19245 

166 

194 

28225 


This table describing the potential mass and cost requirements (as of June 2013) for the 
radioactive isotope components of a hypothetical 300 W RTG is from R. C. Langham. 
“Feasibility study and system architecture of radioisotope thermoelectric generation 
power systems for USMC forward operating bases.,” M.S. thesis, Dept. Sys. Eng., Naval 
Postgraduate School, Monterey, CA, 2013 


11 










































The heat radiated by these materials can be used in a variety of ways to generate 
electrical energy, but, the most direct and thereby one of the more efficient ways, is 
conversion of the blackbody emission of energy proportional to the temperature of the 
source through the use of thermo-photovoltaics. With selective emitters and thermal 
concentration techniques, a material radiating sufficient heat could even be used to 
generate electricity with regular photovoltaics, or a combination of photovoltaics and 
thenno-photovoltaics. 
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III. DIRECT ENERGY CONVERSION 


A. PHOTOVOLTAICS 


Modern photovoltaic cells are essentially large photodiodes intended to maximize 
the photoelectric effect by increasing the area of the cell exposed to incident light and 
produce current through the diode and voltage across the diode’s p-n junction. They are 
generally wide, flat, and thin, with a full surface electrical contact on the bottom—which 
can be mirrored to reflect light back through the cell and generate more current—as well 
as very thin electrodes on the top—to minimize the light blocked. 

As an electrical component, any PV cell has a set of physical characteristics that 
are used to describe and characterize it, and there are generally agreed upon 
characteristics that are used throughout PV research and industry to compare them. The 
more physical of these are the cell’s semiconductor materials, including the impurities 
used to give layers of those materials either an acceptor (P) or donor (N) type, the 
thickness of the P, N, and other layers, as well as the material, thickness, and arrangement 
of the devices electrodes. The electrical characteristics are the cells open-circuit voltage 
(V oc ), short-circuit current (J sc ), IV-curve, fdl-factor, and efficiency. Both the physical 
and electrical characteristics can be used to compare, and through modeling, optimize 
various configurations of PV cells. 

The electrical characteristics of a PV cell are largely based around the cell’s 
efficiency (r|, eta) defined as electrical power out, which is voltage (V) times current (J), 
divided by the power in, which is based on Insolation (I), or wattage per area of light 
incident on the cell, and the area of the cell (A). 


- p 

9 out 


V*J,P = I * A 


( 1 ) 


V oc is the theoretical maximum voltage at which a given PV cell can operate, and 
J sc is the theoretical maximum current at which a given PV cell can operate. 
Unfortunately, neither of these two is a practically useful operating parameter, because a 
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PV cell operating at V oc is also operating with zero current, and a PV cell operating at J sc 
is operating with zero voltage. The practical range of PV cell voltage and current 
operation is described by a graph commonly known as its IV-curve, which details the 
current output from the anode of the cell between zero volts and its V oc , examples of 
which are provided, for several common PV semiconductor materials, in Figure 2. 


Figure 2. IV-Curves for Si, Ge, GaAs, and GalnP Cells. 



This figure describing the basic IV curves for Germanium (Ge), Silicon (Si), Gallium- 
Arsenide (GaAs), and Indium-Gallium-Phosphide (InGaP) and demonstrating the 
variation in J sc and V oc among different semiconductor materials is from P. 
Michalapoulos. “A novel approach for the development and optimization of state-of-the- 
art photovoltaic devices using Silvaco,” M.S. thesis, Depts. of Com. Sci. and Elec, and 
Comp. Eng., Naval Postgraduate School, Monterey, CA, 2002 


Multiplying the operating voltage on an IV-curve by the associated operating 
current results in a graph of power over the voltage range for the PV cell. The maximum 
operating power for the cell, known as Pm 3 x, corresponds to a specific voltage and 
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current, respectively, known as VMax and Jyiax- These are the ideal operating parameters 
for a given solar cell, and Pm 3 x is generally the described P ou t from which the PV cell’s 
efficiency is calculated. A PV cell’s fill-factor is the area of the IV curve accounted for at 
PMax and is suggestive of how ideally a cell operates. The IV curve shown in Figure 3, 
with a graph of Power overlaid on top, graphically demonstrates the relationship among 
these values. 


Figure 3. Sample IV and Power Curves for a Si PV Cell with 
Other Associated Electrical Parameters 


0.030 


0.010 



Current 
♦ Power 


This figure, showing an IV curve overlaid with a Power vs. Voltage curve for a basic Si 
PV cell, as well as the major electrical characteristics considered representative of a cell’s 
overall performance and how those characteristics are calculated, is from data produced 
during early model validation for this research 
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V oc and J sc are somewhat inherent in the semiconductor material or materials 
utilized in constructing the cell. For example, Silicon (Si) cells generally have a Voc of 
0.6 V, and Gallium-Arsenide (GaAs) cells generally have a Voc of 0.9 V. However, even 
these factors can be varied, within the limitations of the material, by changing the 
thickness or doping ratio of the semiconductor material. So, the ultimate design of a PV 
cell can be varied considerably to optimize it for a given application. PV research has led 
to several general techniques employed to optimize a cell’s overall efficiency, largely 
having to do with countering the major sources of efficiency loss for a PV cell. 

One of the greatest potential sources of loss, up to 36%, is reflection of light by 
the surface of the cell [10]. This is typically countered by the use of an Anti-Reflective 
Coating (ARC), which is a thin layer of non-conductive material that is transparent at and 
around the ideal wavelength for that cell. For example, Si02 is one potential ARC for 
silicon solar cells. Though ARCs have been employed in various optics associated 
industries for quite some time, they continue to be a developing area of research in 
association with PVs. Modem ARCs can reduce the light reflected by a PV cell’s surface 
to as little as 5% [10]. One part of the losses associated with reflected light is related to 
light that is at a low angle of incidence to the surface of the PV cell. One solution is the 
use of a textured surface, to alter average reflective angles of the cells surface. Various 
implementations of this have led to PV cell efficiency increases as high as 8% [10]. 

Another major source of loss has to do with the band gap (E g ) of a given 
semiconductor material. The band gap of a material is the energy necessary to sufficiently 
excite an electron in the outer shell of its atomic structure such that that electron is free to 
move through the material. The band gap for conductors is effectively zero and the band 
gap for insulators tends to be much higher than thennal energy or optical sources can 
provide. Semiconductors have a band gap that is relatively easily overcome by thermal 
energy or energy provided by photons; however, there is a minimum energy, the band 
gap (Eg), required to excite an electron for conduction. Si and GaAs have a E g of 1.12 eV 
and 1.42 eV, respectively, which means that a photon with energy lower than that will not 
have sufficient energy to excite an electron into the material’s conductance band. Due to 
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the differences among their band gaps, different semiconductor materials are only 
responsive across a certain range of the EM spectrum above the cutoff threshold set by 
their E g . 

The manufacturing technique addressed to overcoming the limited spectral 
response of various semiconductor materials is the construction of PV cells with several 
layers of materials with different band gaps. These cells are known as dual or multi¬ 
junction PV cells if the layers are physically integrated or as tandem cells if the layers are 
physically or electrically isolated. Figures 4 and 5 show the general construction of a 
multi-junction PV cell, and its associated spectral response. It can be seen that these cells 
have a much wider range of response to EM energy, and are thus more efficient at 
converting incident EM energy over a wide range of wavelengths, such as solar radiation. 
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Figure 4. Sample Multi-Junction PV Cell Design. 



This figure showing the layering employed in a typical multi-junction PV cell, where 
there are three main P-N junctions, the “top,” “middle,” and “bottom” cells, and each is 
separated by a tunnel junction that facilitates the transfer of electric current from one 
junction to another is from P. Michalapoulos. “A novel approach for the development and 
optimization of state-of-the-art photovoltaic devices using Silvaco,” M.S. thesis, Depts. 
of Com. Sci. and Elec, and Comp. Eng., Naval Postgraduate School, Monterey, CA, 2002 
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Figure 5. Spectral Response of a Multi-Junction PV cell, 
Distinguished by Layer 



This figure describing the spectral response range for several different semiconductor 
materials that could be employed in a multi-junction cell, as well as how they overlap and 
could be employed to manufacture a cell with a wider effective band gap (E g ) is from P. 
Michalapoulos. “A novel approach for the development and optimization of state-of-the- 
art photovoltaic devices using Silvaco,” M.S. thesis, Depts. of Com. Sci. and Elec, and 
Comp. Eng., Naval Postgraduate School, Monterey, CA, 2002 


More basic manufacturing techniques to improve efficiency include a full surface, 
mirrored back electrical contact, and the inclusion of a back-surface field. The full- 
surface mirrored back contact allows for low conduction resistance for potential at the 
rear of the cell as well as reflectance of incident light that passes through the cell, such 
that it passes back through the cell and has additional opportunity to excite an electron. 
The back-surface field is a region of higher doping at the back of the cell, meaning there 
are more injected charge carriers in that layer, which helps to increase the cell’s overall 
electrical potential. This layer is in addition to the normal P and N layers of the cell, and 
is typically of the same doping type as the cells lower layer, but with a higher 
concentration of impurities (e.g., P+ if it is an N on P cell or N+ if it is a P on N cell). 
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While the overall efficiency improvement from these techniques is not as high as those 
previously mentioned (1-3% in the case of the mirrored rear contact) [10], it is typically 
sufficiently substantial to warrant the additional associated cost. 

Today there is a significant range of photovoltaic cells in production with an 
equally significant range of efficiencies. Simple cells that use a single, non-crystalline, 
semiconductor material with a single p-n junction, can have an efficiency lower than 8%. 
Multi-junction cells, with a range of semiconductor materials tailored to specific bands, 
placed under solar concentrators have been independently verified to have an efficiency 
as high as 46% [11]. Figure 6 is a regularly updated diagram provided by the U.S. 
National Renewable Energy Laboratory (NREL) that lists a range of photovoltaic 
technologies and their development over time, with the most recent and highest 
efficiencies for each technology highlighted at the right of the figure [1]. 

Figure 6. NREL Best Research-Cell Efficiencies 
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This figure describes the overall best PV cell efficiencies achieved for a given vein of PV 
technology in a given year and is from Best research-cell efficiencies. (2015). National Renewable 
Energy Lab. [Online]. Available: http://www.nrel.gov/ncpv/images/efficiency_chart.jpg 
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The most common PV technologies employ silicon (Si) cells, though this is now 
largely due to their early development and maturity of manufacturing technology, which 
is becoming increasingly diverse and widespread. Though less efficient than Si cells 
during their initial development, gallium-arsenide (GaAs) cells quickly became more 
efficient [1], and became the technology of choice for use in PV spacecraft electrical 
power systems. With the development of thin-film technologies and the incorporation of 
impurities which allow for wider band gap reception (discussed further in Chapter IV), 
GaAs-based cell technologies seem poised to dominate the regions of the PV market 
which are focused on overall cell efficiency with fewer concerns for cost [1]. 

Though current designs are still far from the theoretical limit for PV cell 
efficiency, which is as high as 86% when the Shockley-Queisser limit is applied to multi¬ 
junction cells under concentrated solar radiation [12], the advancements in both design 
and production implemented for modern cells have allowed PV systems to become 
considerably more efficient, cost-effective, and widely used. 

B. THERMOELECTRIC GENERATORS 

Thermoelectric generators (TEGs), are a predecessor to thermo-photovoltaics in 
various direct heat to electricity conversion systems. TEGs will be discussed here to 
emphasize the differences in how they are used in a direct energy conversion system vice 
how thermo-photovoltaics are used. They will not be included as a component of 
modeling or analysis through the remainder of the thesis. 

TEGs were one of the first devices developed for direct conversion of heat to 
electricity, and are still used today. They utilize the thermo-electric effect—which 
includes the Seebeck effect (discovered by Thomas Johann Seebeck in 1821), the Peltier 
effect (discovered by Jean Charles Athanase Peltier in 1834), and Thomson effect 
(discovered by William Thomson, also known as Lord Kelvin, in 1851) [13]. 

The Thomson effect describes how a gradient in current density caused by the 
dissimilarity in temperature at one end of a conducting material and the other will cause 
electricity to flow from one end of the material to the other. The Seebeck effect describes 
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how different conducting materials have different current densities at the same 
temperature, which can result in a voltage flowing between them when they are 
connected [13]. The proportionality of this difference is related by a material’s Seebeck 
coefficient. 

When strands of conducting materials with different Seebeck coefficients are 
electrically connected in an alternating zig-zag pattern, with one side of the zig-zag, 
where the dissimilar conductors meet, thermally connected to a cold plate and the other 
side thermally connected to a hot plate, as seen in Figure 7, a diagram provided by the 
Naval Power Unit in Port Hueneme, CA [14], this will result in a voltage through the 
device. Modern thermocouples often use P and N doped semiconductor materials, which 
can have very high Seebeck coefficients, but very low thennal conductivity, which helps 
insulate the hot plate from the cold plate through the device. 
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Figure 7. Basic TEG design 



THERMAL CONDUCTOR, ELECTRICAL INSULATOR 
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(p) r-TTRE THERMOELEMENT 


This figure describes a basic TEG design and is from “Radioisotope thermoelectric 
generators of the U.S. Navy,” USN, Naval Power Unit, Port Hueneme, CA. Jul. 1978. 


Unfortunately, TEGs are generally not very efficient, particularly with a small 
temperature difference (AT) across the device. Though theoretical efficiency for a given 
device can be complex to calculate, experimental results for a modern commercial TEG, 
obtained by O’Halloran and Rodriguez [15], showed a maximum efficiency of only 
2.22% with an average AT = 68.1 °C. They later stated that other research has shown 
TEGs can potentially achieve an efficiency of 10% with a AT = 500 °C, but even this 
efficiency is a fair bit lower than modern PVs and would require remarkable engineering 
to maintain the required temperature difference. 

Despite their lack of efficiency, TEGs are productive and reliable devices that 
have been used in radioisotope thermoelectric generators (RTGs) for some time. An RTG 
uses the heat of a decaying radio-active element to supply heat to a device that converts 
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that heat to electricity, as seen in Figure 8, normally either TEGs or TP Vs. These devices 
have been employed for electric power production on missions that require an extended 
and reliable energy source in places where PVs cannot practically be employed, such as 
deep space and the deep sea [9], [14]. However, because TPVs are thinner and lighter 
than TEGs and do not require an additional system to maintain the temperature gradient 
across the device—which is not to say they do not require cooling, just that it is less 
complex and intensive—the use of TPVs in an RTG would potentially increases their 
efficiency and reduce their weight and complexity. 


Figure 8. Basic RTG Design. 


Radioisotope Thermoelectric Generator (RTG) 


The RTG converts heat to useful electricity using no 
moving parts. 



Heat source: 
radioactive plutonium 
pellets encased in 
protective materiol 


The heat from the decaying 
plutonium is convened to 
electricity by means of pairs 
of metal plates called 

thermocouples 


This figure shows a modem RTG that has been used in several NASA deep-space 
program power systems; it primarily consists of a radioactive heat source (Pu-238) and 
TEG elements to convert the emitted heat to electricity. The figure is from K. Tate. 
(2011, Nov. 21). Nuclear generators power NASA deep space probes (infographic). 
[Online]. Available: http://www.space.com/13702-nuclear-generators-rtg-power-nasa- 
planetary-probes-infographic.html 
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C. THERMO-PHOTOV OLTAICS 

Thermo-photovoltaics (TPVs) operate in the same manner as PVs, such that all of 
the points and considerations discussed in Section II.A are as applicable to TPVs as PVs; 
however, TPVs are explicitly intended to operate in the infrared, vice visible, wavelength 
range of the EM spectrum. They utilize semiconductor materials that have a lower band 
gap, and are thus responsive to lower energy, longer wavelength photons. 

The infrared range technically includes any wavelength between 0.7 and 
1000 microns, but the practical range for infrared, terrestrially, is from 0.7 to 1.4 microns 
and 1.4 to 3 microns, which are considered the near-infrared (NIR) and short-wavelength 
infrared (SWIR), respectively. At the upper end of this range, at about 2.5 microns, our 
atmosphere is effectively opaque, which means that higher wavelengths are absorbed or 
attenuated by air over a relatively short distance. Thermo-photovoltaic cells do not 
generally operate over that entire range, as they are still limited by the band-gaps of the 
semiconductor materials of which they are constructed. Their intended applications have 
generally been as an internal component of a high heat producing system with a nominal 
temperature considerably higher than 300 K, such as an RTG or propane-gas engine, and 
with a peak wavelength in the IR range itself, near the band gap of the TPV material. 
Research into potential TPV application, such as that by Presby [16] and Davenport [5], 
suggests that employment of TPVs in systems that normally employ TEGs could result in 
considerable efficiency gains. 

To appreciate why TPV cells have not historically been considered viable for use 
in a manner similar to PVs, it is useful to understand the concept of a black-body curve. 
A blackbody curve is the theoretical output, described as spectral energy density at a 
given wavelength throughout the EM spectrum, for a body at a given temperature. For 
example, our sun is generally considered to emit radiation similar to a 5800 K blackbody. 
Example curves for the radiation emitted by a blackbody at several different temperatures 
are given in Figure 9. The curves described are quite pronounced, and have a distinct 
peak wavelength (k Ma x), primarily dependent on the temperature (T) in Kelvin. 
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Figure 9. Blackbody Emission Curves at Various Temperatures. 
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This figure shows the blackbody emission spectra curves at several temperatures over the 
EM wavelength range from 0 to 2 pm, and is from early work in this thesis involving the 
generation of data for a 1300K black body spectrum. 

The peak temperature of a blackbody spectrum can be found using Wien’s 
displacement law, where b is Wien’s displacement constant (2.897E-3 mK): 

4 „, = j ( 2 ) 

The amount of radiation emitted above and below the peak wavelength drops off fairly 
sharply to either side, particularly at higher temperatures. Our sun, again, emitting 
radiation similar to a 5800 K black-body, has a peak wavelength near 500 nm, and the 
majority of the radiation it emits is in the range of 380-750 nm (or 0.3-0.7 microns), also 
known as the range of visible light. This means that the majority of the Sun’s emitted 
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radiation has a much higher energy level than what would be ideal for most TP Vs, which 
have a much lower optimal band gap. Determining the optimal band gap for a PV or TPV 
cell converting energy from a body emitting radiation at a specific temperature is 
relatively simple utilizing Planck’s constant (h = 6.626E-34 J-s) and the speed of light 
(c = 2.998E8 m/s) to find the energy (E) of the photons at the peak wavelength for that 
body; finding the energy, in electron Volts (eV) relative to a wavelength in microns (pm), 
is even simpler: 


E 


he 

T 


;E(eV) 


1.24 

A(jum) 


(3) 


Whether for solar radiation or an enclosed RTG-like device, the optimal band gap 
for a semiconductor material to be used in a PV or TPV cell is one that is just below the 
peak wavelength for the emitting body. Longer wavelength photons in light incident on 
the cell will not have sufficient energy to excite an electron above the band gap of the 
material, and will either be lost or absorbed as heat, while shorter wavelength photons 
incident on the cell will have more than enough energy to excite an electron over the 
band gap of the material, but the additional energy will be lost as heat as well. However, 
there are still significant volumes of energy emitted by a blackbody-like radiator outside 
of the range used to determine an ideal band gap, and use or integration of PVs with 
TPVs could allow for harnessing that energy to improve overall PV cell efficiency. 

Because TPV cells generally operate at a lower band gap, they typically output at 
a lower voltage but higher current. Aside from simple circuit design, placing cells in 
series and parallel, various techniques have been considered for optimization of TPV 
cells. One of these is the subdivision of a single cell into numerous close rows of smaller, 
narrower cells that are connected in series, below the level of the cell as a whole, in order 
to amplify the voltage output of the cell, though at the expense of the cell’s current 
output. Another method has been the use of selective emitters, as suggested by Nefzaoui 
et al. [17], which are materials that emit radiation at specific wavelengths when heated or 
act as filters, only allowing the passage of specific wavelengths, while reflecting others 
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back to the emissive source. In theory, these materials could be designed such that they 
emit radiation specifically at the band gap of the TPV cell material in use. 
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IV. MODELING SOFTWARE AND ASSUMPTIONS 


A. SILVACO ATLAS 

SILVACO has been producing modeling software for semiconductor devices for 
over 30 years [18], and their ATLAS software has been used to model solar cells for 
quite a bit of that time. Its practicality of use for modeling potential solar cell 
technologies was evaluated independently at the Naval Postgraduate School (NPS) first 
by Michalapoulos [19] in 2002, later by Davenport [5], Presby [16], Canfield [20], and 
Utsler [21], It is still in use for various modeling projects, both those relating to PV 
research and, more generally, to just about any electronic circuit component employing 
semiconductor materials. A selection of features and benefits purported by Silvaco for 
use of ATLAS, that are related to this effort, include [22]: 

(1) Electrical, thermal and optical characterization of advanced semiconductor 
devices allows for device performance optimization 

(2) Largest selection of silicon, III-V, II-VI, IV-IV, or polymer/organic 
technologies, including CMOS, bipolar, high voltage power device, 
VCSEL, TFT, optoelectronic, LASER, LED, CCD, sensor, fuse, NVM, 
ferro-electric, SOI, Fin-FET, HEMT, and HBT 

(3) Parallel processing supported on multi-core, multiple processor SMP 
machines and distributed computing 

(4) Addresses challenges of current technologies to help users reduce product 
development time 

(5) Exploration of novel device technologies for next-generation devices 


29 



1 . 


Overall Software Functionality 


Silvaco ATLAS is one component of the Silvaco TCAD suite, and is used for the 
design and operating simulation of two or three dimensional models of semiconductor 
devices. The software employs a system similar to a finite element model to represent the 
physical properties of a device and to propagate the effects of various boundary 
conditions set for the device during a simulation [23]. The majority of initial physical 
characteristics are set by the design of the device and the properties of materials specified 
in the device model that are either drawn from databases built into the software, 
databases that are closely linked to the software—such as optical parameters drawn from 
the Sopra database—or that are specified by the user within the model itself [23]. 

When simulating the operation of a PV or TPV device, we are primarily 
concerned with how the modeling software propagates charge carrier concentration— 
specifically—electron and electron-hole migration, with the associated electric field 
potential—and how it propagates the absorption or transmission of EM radiation—to 
detennine the potential photogeneration for a given device design. For electron and 
electron hole migration, Atlas employs the following equations, respectively, where n or 
p is the carrier concentration at any given node, q is the charge of the carrier, >4J is the 
intrinsic potential, 4> is the potential corresponding to the Fermi level of the material, k is 
Boltzmann’s constant, and T is the temperature in kelvin [23]. 


n = me exp 


kTi. 


( 4 ) 


p = me exp 


-q(y-(/>p) 

m 


( 5 ) 


For absorption or transmission of EM radiation, ATLAS employs a sub-module 
called Luminous. This module performs ray-tracing and calculates the depth of 
propagation of EM waves for wavelengths tested in a given simulation based on the 
material’s refractive index (n) and extinction coefficient (k). Though there are a variety of 
models that can be employed via the Luminous module, simple ray-tracing calculations 
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use the following formula to calculate reflection and transmission, where 0; represents 
the angle of incidence, 0 r represents the angle of reflection, and 0 t represents the angle of 
transmission [23]. 


Or = Oi 

(6) 

m sin Oi = U 2 sin Or 

(7) 


Luminous calculates EM radiation absorption and photogeneration based on the 
following equations, where P represents the ray intensity factor, qo is the material’s 
internal quantum efficiency, y is relative distance along the traced ray, h is Planck’s 
constant, X is the wavelength of the radiation, c is the speed of light, and a is the 
calculated absorption coefficient [23]. 


PA 

= 770 - ae y 

he 

( 8 ) 

Atz , 

a =- k 

A 

( 9 ) 


These equations for charge carrier density, EM wave propagation, and 
photogeneration are examples of the calculations perfonned during an ATLAS simulation 
run to determine the state of charge carrier generation and transport throughout a device, 
as well as overarching electrical properties relating to boundary conditions set for a 
device, such as the voltage at a specific electrode. ATLAS then employs linear and non¬ 
linear solver routines to approximate these conditions throughout a device, and determine 
the properties which are being solved for in any given simulation, such as voltage or 
current output at a specified electrode. 

2. Data Structures and Output 

There are two primary data outputs from an ATLAS model run. The first is the 
actual semiconductor device model, called a structure (a .str file). This is essentially the 
finite element model of the device, with a two- or three-dimensional array of nodes 
representing the devices physical characteristics. The density of these nodes can be varied 
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based on the way the mesh—which is the network of nodes that overlays the material 
layers of the device—is specified, meaning that finer model resolution can be employed 
in areas of particular consideration for a given device. For example, the meshes specified 
for the PV, TPV, and tandem cell models have a higher density of nodes in the emitter 
layers, and at the P-N junctions, but a lower density of nodes in the areas of bulk 
substrate. The second main data file outputs from an ATLAS model run are log files 
(.log); these files store the solutions obtained by the ATLAS solver, which represent the 
resulting optical, electrical, or other properties calculated by the model under set 
conditions. For example, a light beam representing the AMO spectrum, over the center of 
a PV cell’s specified width (a variable in this statement), at a 90° angle of incidence, 
which can undergo reflection by elements within the electronic device, can be specified 
in the solver using a statement similar to the below, where the first term specifies the 
characteristic of the model being described—that this is a “beam”— the second term is an 
ordinal index for that characteristic (for reference by other statements and in case more 
than one of that element type is specified), and the remainder of the terms are specific to 
that element, as previously suggested. 

beam num=l x.origin=$cell_width/2 y.origin=$l_top-10 
angle=90 AMO back.refl 

The solver can then be run with a statement specifying the use of this beam 
model, and telling ATLAS that you would like it to solve for the current generated in a 
PV cell by specific wavelengths of light in the AMO spectrum, in order to generate a 
curve for the cell’s spectral response. Numerous other characteristics can be specified in a 
similar manner, and can be incorporated as a component of the solver’s calculated 
solutions. While a .str file will contain the entire model of a semiconductor device, the 
log file will contain only what is solved for and the relevant conditions for the cell, such 
as an array of values for generated current based on the incident wavelength of light. 

Both types of data file can be viewed in the suite’s companion application 
tonyplot, which will interpret either a .str file into a display of the device, relevant to the 
generated data you wish to view such as layers of different materials, a wireframe 
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representation of the model nodes, or calculated photo-generation under the conditions 
specified during the solver run. It will also generate x-y plots of the data in a .log file for 
analysis of device properties such as spectral response and current vs. voltage (IV) 
curves. Examples of each are provided in the Results section. While no other application 
would likely be viable for tonyplot’s .str file plotting, its somewhat unwieldy interface 
means that its nonnal 2D data plotting capabilities leave much to be desired; however, 
tonyplot also functions for converting the Atlas .log files into normal data file fonnats 
such as comma-separated-value (.csv), which can be used in applications such as 
MATLAB and Microsoft Excel. 

3. Using Deckbuild 

The most widely used operating system, of those on which the ATLAS suite is 
designed to run, are the Microsoft Windows operating systems. Ironically, the normal 
Graphical User Interface (GUI) for semiconductor device modeling, designed to function 
with this software suite, does not work on Microsoft Windows. This is a kn own issue, 
and one which Silvaco seems to have no intention of changing. Thus, the normal 
interface for device modeling on a Microsoft operating system is a script-based interface 
called Deckbuild. While a script-based interface would seem far less intuitive than a GUI, 
the operation of the GUI itself is somewhat redundant as it is essentially used to generate 
the same commands that would be used in the script, in the background, in order to create 
the same device model you would with the script-based interface. While the lack of a 
GUI may seem frustrating at first, and may be particularly challenging for someone with 
little or no programming experience, it eventually allows for a no less, and potentially 
more, consistent model design and testing results, particularly when the user comes to 
understand both the base commands and that the remainder of the ATLAS application 
suite can be integrated into the analysis process through the Deckbuild script interface. 
This method can also better lend itself to integration with an overarching command-line 
interface that could be used to run multiple iterations of a script with varying parameters 
and thus varying models of the device based on those parameters. While this method was 
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not employed in this research, it has been in similar research [21] and has been found to 
be advantageous for the collection of highly granular data. 

The majority of script commands and explanations are outlined in the Atlas User 
Manual [23], along with a great deal of detail on the different models employed by the 
solver, and their related calculations. While a thorough overview or examination of these 
is beyond the scope of this thesis, it is useful to have them readily available for the self- 
education of someone using the Atlas suite, such that they can correctly apply the models, 
solver instructions, and commands relevant to the semiconductor device modeling efforts 
they are conducting. However, it must be understood that—though there is likely overlap 
with other Silvaco software products—the numerous script commands and syntax 
variations used, while logical, are specific to the Atlas suite. Even with a basic 
background in electronics, semiconductor physics, programming, other similarly related 
fields, and even with examples provided with the software, these commands can take a 
great deal of time and practice to learn to employ correctly. For this reason, it may not be 
worthwhile to employ these modeling tools for a more limited effort, without a long-tenn 
gain for the user (e.g., a future career which heavily employs semiconductor device 
modeling), or without significant support by personnel more experienced with the 
software suite. 

B. GAAS PV CELL 

Gallium Arsenide (GaAs) PV cells are historically some of the most efficient 
available, and are commonly used in spacecraft applications. Though even single junction 
GaAs cells are considerably more efficient than Silicon (Si) cells—with even single 
junction single crystalline cells showing an efficiency of about 16%, as seen in Figure 
6—GaAs is also a common base or substrate in multi-junction cells. One of the reasons 
for this is that GaAs, particularly in combination with Indium (In), Phosphorous (P) and 
some other impurities, has a tunable band gap. 

While the potential for multilayer InGaAsP devices has certainly been considered 
and is being explored in modern research and design [21], the intent of this effort was to 
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determine the simple potential efficiency improvement of the incorporation of a PV cell 
with a TPV cell, which will be best achieved by eliminating extraneous variables such as 
variability of impurities, which could result in band gap variability as well. As a 
common, well-modeled PV semiconductor material, with a sharp spectral response 
termination at 0.9 pm, a simple, single-junction GaAs cell, with distinct N and P doping 
regions, should serve well for this purpose. A sample GaAs PV cell configuration is 
shown in Figure 10. 


Figure 10. Sample GaAs PV Cell Plot 
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This plot describes the materials used in a sample tandem PV cell design, as well as the 
thickness (pm) and relative positioning for layers of each material. These thicknesses are 
not representative of any cell modeled used during this research; they are exaggerated for 
the sake of visibility, as normal layer proportions make it difficult to separately 
distinguish layers in a static image. 


To this end, research was conducted to determine common parameters for a GaAs 
PV cell. Even with this restricted scope, it was found that there is still considerable 
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variability in the parameters employed in manufacturing a simple GaAs PV cell. The two 
most commonly varied parameters are substrate thickness, and doping profile. Having an 
expectation that the overall thickness of the GaAs top cell, in a tandem style 
configuration, would have a dominant effect on the efficiency of the lower cell, and 
knowing that the substrate layer of a PV cell is normally the thickest, it was decided to 
vary the substrate thickness of the cell across a relatively broad range (from 5-500 pm), 
to determine its effect on overall efficiency. The effects of doping profile can be equally 
pronounced, though to limit the problem space it was decided that two basic categories of 
doping profile, with unvaried constituent values, would be used for comparison: a low 
doping profile, and a high-doping profile, the values of which are provided in Table 3. 
Basic modeling efforts, prior to commencement of final experimental modeling, 
demonstrated that the inclusion of an ARC layer and a BSF layer produced model results 
more consistent with realistic values, so these were employed in final experimental 
modeling. Though both N on P and P on N methods have been used for GaAs PV cells, it 
was chosen to employ an N on P design with a P+ BSF, because this configuration 
appears to be more common. 


Table 3. Doping Profiles for GaAs PV Cell. 



Low Doping profile (LD) 

High Doping profile (HD) 

Top (N) Layer 

1.00 E17 

1.50 E18 

Substrate (P) Layer 

1.00 E16 

7.00 E17 

BSF (P+) Layer 

1.00 E17 

8.00 E18 


C. GASB TPV CELL 

Contrary to the decision supporting the use of GaAs for the PV cell modeling, it 
was decided to use GaSb for TPV cell modeling because its coverage in research has not 
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been as thorough, but because it should serve equally well as an example of a TPV cell 
material as GaAs serves as an example of a PV cell material. 

Because the majority of GaSb TPV cells found in research utilized a P on N 
design, this was chosen, with an N+ BSF, as the general construction for the TPV cell. A 
sample plot of a GaSb TPV cell is shown in Figure 11. A configuration of contacts and 
ARC similar to the PV cell were employed for testing of the single TPV cell. 


Figure 11. Sample GaSb TPV Cell Plot (with Mesh Model Overlay) 



Cell Width (pm) 


This plot describes the materials used in a sample tandem PV cell design, as well as the 
thickness (pm) and relative positioning for layers of each material. It also shows an 
overlay of the mesh used to model the nodes of the cell, as well as the highlighted anode 
and cathode contacts and a purple line representing the p-n junction. Again, the 
thicknesses shown are exaggerated for visibility and are not representative of any cell 
model used in this research. 


Similar to the PV cell, two doping profiles were chosen for the TPV cell, one 
found experimentally to be sufficiently effective, relative to expected TPV efficiency 
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values—the Low Doping profile—and the other extracted from Davenport’s research 
detailing the optimization of a GaSb TPV cell for a 1300 K blackbody [5]—the High 
Doping profile. 


Table 4. Doping Profiles for GaSb TPV Cell. 



Low Doping profile (LD) 

High Doping profile (HD) 

Top (P) Layer 

1.00 E17 

9.80 E19 

Substrate (N) Layer 

3.50E17 

1.00 E20 

BSF (N+) Layer 

1.00E18 

1.00 E20 


D. TANDEM GAAS-GASB PV-TPV CELL 

The basic design for the tandem GaAs-GaSb PV-TPV cell is simply the PV cell, 
on top of the TPV cell, with the two electrically isolated from each other by a layer of 
“gap” material. A simple example model is shown in Figure 12, with the thickness of the 
layers exaggerated for visibility; in most of the modeling runs, the thickness of one layer 
relative to the others—as well as the contacts, ARC, and gap material—makes it difficult 
to visibly distinguish them from each other in a static image. Though tandem designs are 
not popularly employed in modem PV cell manufacturing, they were employed during 
the early years of PV cell development and were able to produce devices of considerably 
higher efficiency than other PV cells of the time in a manner similar to that of modern 
multi-junction PV cells, where the upper layer is intended to absorb photons at shorter 
wavelengths, successive lower layers are intended to absorb photons with longer 
wavelengths, and the material used at each layer possesses a band gap which best fits to 
the radiation which penetrates to its depth in the cell. The overall complexity and cost of 
manufacturing a tandem cell, along with the discovery of multi-junction implementation, 
was likely the driving factor in the reduction of interest in the design; however, modeling 
and simulation software allows for continued exploration of this design, with relatively 
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few concerns about cost, and the complexity of a modelled tandem design is arguably 
less than that of a multi-junction design, particularly in regard to degree of computation 
required by the model. 


Figure 12. Sample Tandem PV-TPV Cell Plot 
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This plot describes the materials used in a sample tandem PV-TPV cell design, as well as 
the thickness and relative positioning for layers of each material. Once more, the 
thicknesses shown are exaggerated for visibility and are not representative of any cell 
model used in this research. 


Considerations were made relative to the ordering of the N and P layers for the 
tandem cell overall, where, in a multi-junction cell, the ordering is required to be 
consistent to allow for the implementation of a tunnel junction between the overall cell 
junction layers, to allow the passing of current from one cell junction through to the 
other. However, research has shown the similarity of function between N on P and P on 
N implementations and both the electrical isolation of the cells from each other, in the 
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tandem model, and the ability to calculate potential overall values for a tandem or multi¬ 
junction cell, given these results, proved this to be an unnecessary consideration, which 
could have doubled the potential problem space for the tandem cell. 

Other considerations for the design of the tandem cell model included the 
arrangement of electrical contacts between the cells, the thickness and type of gap 
material used, and the use of ARC between the cells. Given that upper electrical contacts 
are a major source of loss in an individual cell, due to the reflection of incident EM 
energy, it was decided that the contacts would be implemented at an appropriate 
thickness and width for a normal PV cell electrical contact, but that they would be 
positioned directly above and below each other to minimize this source of loss. 
Additionally, while a normal PV cell has thin, finger-like upper contacts to minimize 
reflected light, but a full surface back surface contact to maximize potential photon 
reflection back through the cell, it would not have been practical to employ a full surface 
contact on the upper of the tandem cells, as it would have blocked incident IR radiation 
from penetrating to the lower cell. Thus both the upper and lower contacts on the PV cell, 
as well as the upper contact on the TPV cell, are of a finger design, and only the bottom 
contact of the TPV cell is full surface. 

The gap material used was TiCE, which is both transparent and non-conductive. 
This allowed for the passage of IR and any remaining visible light energy through the top 
cell to the bottom, while preventing any potential passage of electrical current. In the 
majority of testing, the thickness of the gap layer was held at 1 pm, to maintain 
consistency of results relative to variation of the substrate thickness of the upper cell. 
However, in later tests of thinner upper cells, detailed in Chapter V, the gap layer 
thickness was significantly reduced to help eliminate its considerations as a potential 
source of energy loss between the upper and lower cells. Finally, ARC was employed 
between the cells in the majority of modeling runs, with its thickness adjusted to tune 
absorption for longer wavelengths of EM radiation for the lower cell. In the thin cell 
modeling runs, similar to the gap material, its thickness was adjusted on the upper cell 
and eliminated on the lower cell, to reduce its consideration as a potential source of loss. 
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Efficiency (%) 


V. RESULTS 


A. GAAS UNDER AM1.5 AND AMO 

Both LD and HD cell models were run over 5-500 |am substrate thickness to 
determine overall expected performance for a GaAs cell under AM 1.5 and AMO 
spectrums and compare to real-world performance values. Full electrical performance 
data for each run is detailed in Appendix A, and Figure 13 shows a comparison of results 
based on overall cell efficiency relative to substrate thickness. 


Figure 13. Comparison of Low Doping and High Doping GaAs PV Cell 
Efficiencies Under AM 1.5 and AMO Spectrums 



I LD: AMI.5 
I - LD: AMO 
►--- HD: AM1.5 
> HD: AMO 


Highest perfonnance, 14.24%, was shown by the LD profile cell under AM1.5 at 
a substrate thickness of 35 |im. Of interest is that both LD and HD profiles performed 
less efficiently under AMO than under AM1.5—by approximately 1.5%—which is likely 
attributable to the bands of EM energy lost to the atmosphere under AM 1.5 being outside 
of the spectral response range of this GaAs cell. Unfortunately, for spacecraft power 
systems our primary interest is efficiency under the AMO spectrum, for which the highest 
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result, 12.52%, was again shown by the LD profile at a thickness of 35 |am. A plot of the 
cell model, the cell’s spectral response curve, and the cell’s IV and power curves are 
shown in Figures 14, 15, and 16, respectively. Though both cells showed respectable 
performance for a single junction mono-crystalline PV cell, under both AMI.5 and AMO, 
they did not show performance comparable to modern single-junction GaAs cells— 
approximately 18-29% [11]. It is believed that the simplicity of the PV cell model design 
is the likely source of this difference in performance. While it was intended to down- 
select to the higher performing doping profde for use in modeling a tandem design, the 
overall similarity among efficiencies for each profile, and a sufficient allowance of time, 
led to the testing of both doping profiles in the tandem design. 

Figure 14. AMO Optimal GaAs PV Cell Model 
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Figure 15. AMO Optimal GaAs PV Cell Spectral Response 
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Figure 16. AMO Optimal GaAs PV Cell IV and Power Curves 
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B. 


GASB UNDER AM1.5 AND AMO 


Both LD and HD cell models were run over 5-500 pm substrate thic kn ess to 
determine overall expected performance for a GaSb cell under AM 1.5 and AMO 
spectrums and compare to real-world perfonnance values. Full electrical performance 
data for each run is detailed in Appendix A, and Figure 17 shows a comparison of results 
based on overall cell efficiency relative to substrate thic kn ess. 


Figure 17. Comparison of Low Doping and High Doping GaSb PV Cell 
Efficiencies Under AM 1.5 and AMO Spectrums 



LD: AM1.5 
— * — LD:AM0 
HD: AMI.5 
o HD: AMO 


The GaSb TPV cell design showed increasing efficiency with substrate thickness, 
well past the 500 pm range. Though this makes it difficult to detennine overall highest 
performance, based on substrate thickness, it can be seen that the HD profile, overall, 
performed considerably better than the LD profile—by 6.5% and 7.5% under AM1.5 and 
AMO respectively. Due to this drastic difference in efficiency, and for the opportunity to 
restrict the problem space to a reasonable size, the LD profile was down-selected for use 
in the tandem design, in favor of the HD profile. 

It can also be seen that there is little variation in efficiency based on substrate 
thickness within an individual doping profile under a given spectrum, the greatest being 
within the HD profile under AMO, where a 5 pm substrate thickness shows an efficiency 
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of 17.82% and a 500 |im substrate thickness shows an efficiency of 18.16%—a 
difference of only 0.34%. For this reason, it was chosen to fix the substrate thickness of 
the TPV cell in the tandem design to a thickness somewhere near the “knee” of the curve, 
at the lowest thickness where the increase in efficiency, relative to substrate thickness, 
begins to level off. Though somewhat arbitrary, this, again, allowed for some restriction 
of the problem space for the tandem design; the chosen thickness was at 50 |am, which 
showed a single-cell efficiency of 18.04%. The cell model, spectral response, and IV and 
power curves for this cell are shown in Figures 18, 19, and 20, respectively. 

Figure 18. AMO Optimal GaSb TPV Cell Model 
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Figure 19. AMO Optimal GaSb TPV Cell Spectral Response 
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Figure 20. AMO Optimal GaSb TPV Cell IV and Power Curves 
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Another important consideration is the relative performance of the GaSb cell to 
the GaAs cell. Contrary to the decreased performance of the GaAs cell between AMI.5 
and AMO, the GaSb cell showed an increase in performance between AM 1.5 and AMO of 
about 2.5%. Again, this can likely be attributed to the differences in spectral distribution 
due to atmospheric absorption, but, in this case, those areas of the EM spectrum absorbed 
by the atmosphere are within the spectral response range of GaSb, where they were not 
with GaAs. Also, though the LD profile, at no point, performed better than the GaAs cell, 
under either spectrum, the HD profile performed better under both AM 1.5 and AMO than 
the peak perfonnance of the GaAs cell. Compared to the 12.52% efficient GaAs 
performance under AMO, the chosen thickness of HD GaSb, at 50 |am, showed an 
efficiency increase of 5.53%—an improvement of 44%. Tough it was important to 
examine the modeled performance of single-cell GaSb, for comparison to its performance 
in a tandem design, because GaSb cells have not commonly been considered as a PV 
material, by themselves, it is difficult to relate these results to anything found in the 
literature. 

C. GAAS-GASB TANDEM CELL UNDER AMO 

Once the expected perfonnance of the individual cell models was established, the 
models were combined to determine a single, tandem cell, with the GaAs PV cell on top, 
the GaSb TPV cell on the bottom, and a layer of Ti02 gap material between to electrically 
isolate the cells from one another. Again, while the HD profile and a set substrate 
thickness of 50 |im were chosen for the GaSb TPV cell, the substrate thickness of the 
GaAs PV cell was varied to determine overall impact on the efficiency of the TPV cell, 
and in search of a near optimal cell thickness, based on overall PV-TPV efficiency. Tests 
for these cell models were initially run under the AMO spectrum only. 

Performance of both the Upper Cell (GaAs, PV), abbreviated as “U” in Figures 21 
and 22, and the Lower Cell (GaSb, TPV), abbreviated as “L” in Figures 21 and 22, were 
independently recorded over substrate thicknesses originally varying between 5 and 
500 |am, as with the GaAs and GaSb single-cell models. These results were then used to 
mathematically calculate the potential performance for either a Tandem or Dual-Junction 
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design. Modeling runs for upper cell substrate thickness less than 5 microns, considered a 
thin cell design for the purposes of this research, were later run, and the overall efficiency 
results obtained in the same way as for the thicker substrates 

Figure 21 shows the overall results for these calculations when substrate thickness 
is varied from less than 1 pm to 500 pm. Over the entire range of upper cell substrate 
thickness, the energy absorbed by the upper cell significantly reduced the efficiency of 
the lower cell; however, over the majority of the range, the lower cell’s efficiency was 
never less than 1.86%. Given the original GaSb TPV cell efficiency was calculated at 
18.04%, this suggest that approximately 10.3% of the incident electromagnetic radiation, 
a significant portion, penetrates to the depth of the TPV cell and is available for 
conversion to electrical energy. As can be seen in Figure 21, the lower cell efficiency 
values for both LD and HD profiles of the upper cell overlap over the majority of the 
range, suggesting that differences in upper cell doping profile do not have an overall 
significant impact on energy available to the lower cell. The perfonnance of the upper 
cell was similar to the results obtained during the single GaAs PV cell modeling, with the 
highest performance being 13.50% by the LD profile at a substrate thickness of 2 pm. 

Unfortunately, the conjunction of these two disparate power elements resulted in a 
significant overall decrease in efficiency over the majority of the range of upper cell 
substrate thickness—from approximately 5 to 500 pm. With the Tandem results, the 
calculations are based on use of power sources of different voltage and current ratings in 
either a parallel or series circuit, whichever could produce the greater power output. With 
the Dual-Junction results, the calculations are based only on the potential result of output 
from power sources of different voltage or current ratings in a series circuit. In a parallel 
circuit, the overall output is restricted by the lowest voltage rating of the given elements, 
and the current ratings are added. With a series circuit, the overall output is restricted by 
the lowest current rating of the given elements, and the voltage ratings are added. With 
the Tandem results, the overall voltage was significantly restricted by the lower cell, and 
with the Dual-Junction results, the overall current was significantly restricted by the 
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lower cell. Thus in each configuration, the overall Tandem cell was, on average 5.3% less 
efficient than the upper cell efficiency alone, an overall loss of over 30% in most cases. 

However, this trend reverses for upper cell substrate thicknesses of less than 
5 |am, as can be surmised from the significant changes in data on the far left range limits 
of Figure 21, and as is shown in detail on Figure 22, which describes the overall Tandem 
and Dual Junction efficiency results for LD and HD profiles. Overall efficiency results 
for the thinner top cell substrates increase to matching, and even exceeding, the results 
for the upper cell alone. It is expected that as the upper cell thins it becomes effectively 
transparent to greater and greater quantities of the incident EM radiation. Though this 
leads to decreasing efficiency for the upper cell, due largely to decreased current 
production, it means that increasing quantities of this energy pass through to the lower 
cell, which, in turn, results in an increase in efficiency based on increased current 
production. These shifts can both be seen in Figure 22 and can be perused in greater 
detail via the data provided in Appendix A. As the upper cell current decreases and the 
lower cell current increases, they eventually come to a point where the current for both 
cells is matched. This is the optimal configuration for power production in a Dual- 
Junction cell design, and it is at this point in both the LD and HD profiles where we find 
their respective highest efficiency ratings. The higher of the two, 15.08%, is once again 
exhibited by the LD profile and is at a substrate thickness of 0.164 qm. The cell model, 
spectral response curves, and IV and power curve for this cell, based on Dual-Junction 
calculations are shown in figures 23, 24, and 25, respectively. It should be noted that the 
Tandem and Dual-Junction results were seen to strictly overlap, for both LD and HD 
profiles, except in the case of the thinnest substrate thicknesses tested, which were near 
the apparent failure point for the upper cell. 
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Figure 21. Upper Cell (GaAs, PV), Lower Cell (GaSb, TPV), and Overall Tandem 

Design Efficiency for LD and HD Profiles. 
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Figure 22. Tandem and Dual-Junction Efficiency for LD 
and HD Profiles (Thin Upper Cell) 
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Figure 23. Optimal GaAs-GaSb PV-TPV Cell Model 
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Figure 24. Optimal GaAs-GaSb PV-TPV Cell Spectral Response 
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Figure 25. Optimal GaAs-GaSb PV-TPV Cell IV and 
Power Curves (Dual-Junction) 
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Of note in Figure 23 is that the relative thickness of the upper cell layers to the 
lower cell make it very difficult to discern in the cell model. Similar disparity can be seen 
in the multi-junction cell example diagram in Figure 4, where the base and substrate 
layers of the design are considerably thicker than those of the upper, emitter, layers. It is 
expected that, as an example, this figure is actually not to scale, and the disparity would 
be considerably more extreme in a true design, which would be much more similar in 
scale to Figure 23. 

Also of note, in Figure 24, is the overlap in spectral response curves for the upper 
cell—GaAs, listed as anode current—and lower cell—GaSb, listed as emitter current. 
While the spectral response of the upper cell cuts off sharply at about 0.9 pm, as expected 
of GaAs, it also decreases more rapidly approaching this point than is seen in Figure 15. 
Conversely, the lower cell shows very little spectral response in the 0.3 to 0.6 [am range, 
but its spectral response then increases relatively rapidly, as the upper cell spectral 
response decreases, to the point where the upper cell ceases to respond, and the lower cell 
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recovers what it can of the remaining EM radiation. As was suggested earlier, the 
decreased spectral response of the upper cell is likely indicative of the lesser current 
production as the cell’s thickness is insufficient to interact with longer wavelengths of 
radiation, which then simply passes through to the lower cell, resulting in its increased 
spectral response. 

Two changes that occurred between the thin cell designs and the remainder of the 
test range must be noted. First, in the majority of the test range, both the thickness of the 
electrodes and of the gap between the upper and lower cells was held at 1 pm. For the 
thin upper cell tests, both the gap and electrode thicknesses where reduced to 0.1 pm. The 
combined decrease in both gap and electrode thickness equated to a total thickness 
reduction, between the upper and lower cells, of 2.7 pm. Second, for the majority of the 
test range, ARC was applied to both the upper and lower cells, with the upper ARC 
thickness being optimized for X = 0.6 pm and the lower cell ARC being optimized for X = 
1.0 pm. For the thin cell tests, there was no ARC applied to the lower cell, and the 
thickness of ARC for the upper cell was optimized for X = 0.5 pm. Though it is difficult 
to quantify the overall impact of these changes, they were made based on a minor overall 
improvement in cell performance seen during various model test runs. It is expected that 
the overall decrease in gap thickness, relative to the thickness in upper cell substrate, may 
have contributed to the majority of this improvement. 

D. GAAS-GASB TANDEM CELL UNDER 2000 K BLACKBODY 

SPECTRUM 

Though the focus of this research has been examination of potential perfonnance 
improvements for a PV cell under AMO by the integration of TPV technology, the 
conjunction of PV and TPV technology could be seen as having the potential to act as a 
hybrid power element capable of producing power relatively efficiently under AMO or by 
utilizing a lower temperature blackbody spectrum such as would be emitted by an RTG. 
To examine this possibility, the cell design optimized for AMO was simulated under an 
approximate 2000k blackbody spectrum—similar to what would be emitted by an RTG— 
to detennine the potential improvements to a TPV design by the integration of PV 
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technology. This design was then varied in a manner similar to that used for rough 
optimization of the tandem cell under AMO to determine a range of potential performance 
for a GaAs-GaSb PV-TPV tandem cell employed in an RTG. 

As calculated by Wien’s Displacement law, a 2000 K blackbody produces a 
spectrum with a peak wavelength at 1.4485 pm, which is slightly higher in energy than 
the Eg wavelength of GaSb, at 1.7 pm. This is a fairly optimal temperature for the GaSb 
TPV cell in the tandem design. However, it is well below the E g wavelength for GaAs, so 
the amount of energy proportionally available to the GaAs PV cell will be considerably 
less than was available under AMO. This is almost the reverse of the AMO case, where 
the Eg wavelength of GaAs (0.87 pm) was slightly below the dominant wavelength of the 
AMO spectrum (0.499 pm), allowing it to absorb the majority of the incident photonic 
radiation, but still allowing a considerable portion of the longer wavelength radiation to 
pass through for absorption by the GaSb TPV cell. The expected results of this situation 
are that the upper PV cell should be able to absorb whatever is available in the smaller 
wavelength portion of the spectrum, but the majority of the radiation should pass through 
to the lower TPV cell. As can be seen in Figure 26, the modeled spectral response of one 
of these cells shows this is exactly what would occur. 
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Figure 26. GaAs-GaSb PV-TPV Tandem Cell Spectral Response Under 2000 K 

Blackbody Spectrum 



Unfortunately, as the data at the end of Appendix A indicate, and as is shown in 
Figure 27, the amount of current available from the upper cell, and thus its efficiency, 
is far below that for the lower cell, resulting in a considerable overall reduction in 
efficiency in a dual-junction implementation of this design. A tandem implementation 
shows an overall marginal improvement—an average increase of approximately 0.8%— 
over the performance of the lower cell alone. In the AMO case, the current output of the 
upper and lower cells was close enough that minor adjustments in thickness allowed for a 
balance in the performance of the upper and lower cells, based on output of electrical 
current, which led to a considerable overall improvement in efficiency using dual¬ 
junction calculations. In this case, the difference in current output is so large that the 
performance of the lower cell would have to be severely degraded in order to obtain a 
similar match, which would result in a considerable overall loss in efficiency. While the 
tandem design does show improvement, the overall output increase of approximately 3% 
may not be worth the added cost and complexity of the design. 
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Figure 27. Upper Cell (GaAs, PV), Lower Cell (GaSb, TPV) and Overall Tandem 
and Dual-Junction Efficiency Under 2000 K Blackbody Spectrum 
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VI. CONCLUSIONS AND FUTURE WORK 


A. OVERALL PV-TPV PERLORMANCE 

While the overall performance of the thicker upper cell designs was discouraging, 
even though indicative of potential efficiency gains, the overall perfonnance of the thin 
upper cell designs was decidedly encouraging, showing great promise for potential 
efficiency gains. Given the overall highest efficiency of both the single cell GaAs PV 
design and the overall highest efficiency of the GaAs-GaSb PV-TPV design, utilizing 
dual-junction calculations—12.52% and 15.08% respectively—we see an increase in 
efficiency of 2.56%, which equates to a 20.4% increase in power output. 

These results are comparable to those found by Andreev, et al. in their 
examination of solar cells based on GaSb, where actual manufactured cells utilizing 
GaSb were examined “behind a filter based on the GalnP/GaAs heterostructure” of a 
potential upper cell in a tandem design [24], Their results showed a potential efficiency 
for the lower cell to be between 4.5% and 6.5% at solar concentrations of 50 to 750 suns 
depending on the material used for the upper cell and its doping profile [24], The same 
study showed tests of a GaSb cell “without a wide-gap cell” covering it to have 
photocurrent density nearly 2.1 times greater than when covered [24]. The range of lower 
cell efficiency results obtained with thin cell models in this study—between 4.4% and 
6.26%—compares favorably with these results. However, the overall expected efficiency 
and photocurrent density for a single GaSb cell in [24] suggest a single cell efficiency of 
no greater than 13.15%, which is quite different from the maximum GaSb single cell 
efficiency of well over 18% found in this study. Though this similarity in tandem design 
results is, again, promising, the suggested difference in overall GaSb cell efficiency 
warrants further investigation. 

It is important to note that the base efficiency of the GaAs PV cell used in these 
models, while comparable to early GaAs PV cell designs, is not comparable to modern 
GaAs PV cell designs. Modem GaAs cell efficiencies are considerably higher, especially 
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for thin film cells which exhibit cell efficiencies in excess of 28% [11]. The manufacturer 
cited for producing thin-film GaAs cells at these efficiencies is Alta Devices [11]. Alta 
Devices employs metalorganic chemical vapor deposition (MOCVD) to manufacture 
these thin film cells [25]. This technology has the potential to manufacture devices at 
thicknesses approaching those of the thin upper cell model designs near which the 
optimal PV-TPV model efficiency was seen in this research [25]. The data in this 
research suggest that use of a TPV cell in conjunction with a modem, higher efficiency 
PV cell design such as this has the potential to result in a PV-TPV cell with efficiencies 
between 30% and 34%. 

The 2000 K blackbody emitter results require some additional consideration in 
conjunction with the AMO results. While the overall performance change for the modeled 
tandem cell design under this spectrum showed little improvement, it is important to 
consider these results relative to both the band gaps of the materials utilized and the 
chosen emitter spectrum. The AMO spectrum was utilized for the majority of the 
modeling tests because it represents the nonnal spectrum available for PV power 
production in the space environment. The 2000 K blackbody spectrum was chosen 
because it represents a sufficiently close approximation to the spectrum emitted within an 
RTG. Neither of these spectrums was chosen to match the ideal potential power output 
for a PV-TPV design utilizing GaAs and GaSb. This is normally the case with renewable 
energy sources, where there is an existing source of energy, and a power production 
system is designed to draw from that, not the other way around. If it were possible to shift 
this design philosophy, such that a power source was designed to emit EM radiation at a 
wavelength that produced ideal results for a given cell design, then an emitter with a 
wavelength closer to 1.3 pm, which would have considerably greater overlap with the 
spectral response range of GaAs, could produce much greater results. This wavelength 
corresponds to about 2230 K, which is not much higher than the emitter temperature 
considered, thus, with the correct blend of radioactive materials in an RTG, this 
temperature may be possible to achieve. 
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B. DUAL-JUNCTION VS. TANDEM PERFORMANCE 

Calculations were performed, given the separate electrical properties determined 
for the upper and lower cells, which would represent the cell’s performance in either a 
dual-junction or tandem design. Because the design of this cell was not directly in line 
with a true dual-junction cell model, the calculated values cannot be seen to be as realistic 
as those that would be produced by a more true dual-junction PV-TPV model, however, 
they are sufficiently accurate, given the detennined electrical properties of each cell, to 
be representative of the values that could be expected from such a model, and they offer 
the benefit of being able to independently examine the electrical properties of distinct 
layers in a dual or multi-junction design, with the modeling software utilized. 

Contrary to their comparison in real-world PV cell design and manufacturing, 
where the two are quite different, distinguishing between these two potential 
configurations, based on the model utilized, is somewhat arbitrary. In this case, with the 
calculations for a dual-junction design, the two cells are treated as though they must be in 
series (as would be the case with a real-world design involving a tunnel junction), which 
means that the voltage output of the two cells is the sum of that from each individual cell, 
while the current output is restricted to the lesser of the two. With the calculations for a 
tandem design, which this model more closely approximates, the cells could potentially 
be configured in series or in parallel, where the series configuration results would directly 
match those for a dual-junction design, but the parallel results would be based on the 
lower of the two voltages between the cells, but the current output would be the sum of 
the two cells. Thus, while the parallel configuration for a tandem design could potentially 
have optimal results quite different from those of a dual-junction design, the optimal 
results for a series configuration are shown to be nearly identical, based on the results of 
this research. Data listed in Appendix A suggest further similarity among these results, 
for the models utilized, where the efficiency results for the thin cell models tested were 
nearly identical for both tandem and dual-junction calculations. 

Only for the thicker substrate designs, where the lower cell efficiency was 
significantly reduced, are there substantial differences among the tandem and dual- 
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junction calculation results. In these cases, the Tandem designs prove more efficient—by 
between 1.21% for the least difference in results with the HD profile and 3.11% for the 
greatest difference in results for the LD profile. However, the overall lower efficiency of 
these designs is not likely to warrant consideration of these differences. 

C. PV-TPV OPTIMIZATION 

Basic optimization methods were utilized to find a rough peak of performance for 
a dual junction or tandem PV/TPV solar cell based on the developed model. The simple 
method utilized was to determine the cells maximum power output and efficiency over a 
range of values, plot a curve of these values based on top / PV cell thickness to determine 
a relative peak, perfonn a finer grained range of calculations in the vicinity of that peak, 
and repeat until there was little or no variance in the highest value determined. 

Other optimization methods have been used in the past to find a greater or even 
absolute theoretical peak of perfonnance for various research explorations, including for 
the parameters of PV and TPV cells. Future research utilizing techniques such as genetic 
algorithms or simulated annealing, in conjunction with the developed Deckbuild scripts, 
could produce considerably more optimal results than those found during the course of 
this research. In past NPS theses, MATLAB has been used in conjunction with Deckbuild 
to apply similar techniques to explore various optimization problems [19], [21]. The 
flexibility and effectiveness of such methods could allow for variation of a considerably 
greater number of both PV and TPV cell parameters, such as doping and thickness of 
individual layers, and could also be used to determine the variation of optimality across a 
range of environmental characteristics. 

D. TEMPERATURE AND RADIATION EFFECTS ON A TANDEM 

PV-TPV CELL 

The space environment is considerably more extreme in temperature than what 
was represented in the models used for this research (an unvarying temperature of 300 
K), and it is known that temperature variation can have a great deal of impact on PV cell 
performance. Various statistical data taken from real PV cells at varying temperatures has 
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shown a trend that temperature increases of one °C generally correlate to a voltage 
decrease of about 2 mV and an amperage increase of about 0.6 inA, normally resulting in 
an overall slight loss in efficiency with increases in temperature. While space itself is 
often considered to be very cold, any object in space that is consistently exposed to solar 
radiation, particularly those that are dark in color, such as PV panels, heats up very 
quickly before reaching a fairly high equilibrium temperature, normally 100 °C or higher. 
Though the correlated 0.2 V loss for this may not seem overly significant, it is important 
to recall that the V oc for the upper cell in this study was just under 1.0 V—thus suggesting 
a potential voltage reduction of over 20%—and the V oc for the lower cell was just under 
0.5 V—thus suggesting a potential voltage reduction of over 40%. An additional 
complication might be that the shading of the upper cell could reduce the resulting 
temperature increase in the lower cell, which would also be nearer the backside of the 
panel and thus cooler, resulting in lower voltage losses, but that the increased heat and 
lesser though still increased temperature would result in substantial gains in current 
production. The temperature impacts to a tandem cell design in a space environment are 
potentially quite complex, and warrant careful consideration when attempting to optimize 
such a cell for this environment. 

Also, the degree of radiation present in the space environment is considerably 
higher and more complex than was represented in the models used for this research— 
which is to say, none. As was expressed in Section II.B., high energy particles and 
photons can have relatively extreme negative impacts on electrical or electronic 
components. While the simplicity of a PV cell, relative to other electronic components, 
generally minimizes the appearance of these effects, they do still occur with sufficient 
frequency to be of concern, and the related overall degradation of cell performance due to 
these effects is considerably more evident in the space environment than it is terrestrially. 
Exploration of the effects of various types of radiation on a PV-TPV tandem cell would 
allow for better understanding of the chance of these impacts causing a critical device 
failure, and would also provide some appreciation for the expected degradation of cell 
performance over time, one of the primary considerations for designing the power budget 
for a long-life system such as most spacecraft. 
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Though it is likely not the only application to do so, SILVACO Atlas does 
provide the capability to explore either or both of these concerns, and the scripts 
developed during this research could be employed for this purpose with fairly minimal 
adaptation. 

E. TEG INTEGRATION 

One of the early concepts behind this thesis was the potential integration of 
various types of direct energy conversion devices, to include thenno-electric generators. 
While, by themselves, they are fairly inefficient, without a high temperature gradient, 
both previous research and this thesis have shown that the proper balance of various 
energy generation techniques can potentially lead to a similar or greater overall efficiency 
than use of a single technique alone. 

Though time did not allow for investigation of this type, it is possible that the 
integration of a TEG layer, in a similar manner to the TPV layer integrated with the PV 
layer in this thesis, could result in additional energy production, given that the “hot foot” 
of the TEG is a PV/TPV array facing the sun and the “cold foot” of the TEG is the 
shadowed side of the array directed at the vacuum of space—normally used for 
spacecraft cooling—or an alternative, even terrestrial, cooling system. 
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APPENDIX A. DATA 


SILVACO ATLAS Modelling Stats for GaAs solar cell (N on P on BSF-P+)(AM1.5) 
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l.OOE+17 

Si02 

0.6000 

1.9672E-02 

8.5955E-01 

1.4O06E-O2 

7.5951E-01 

1.8441E-02 

82.83% 

14.01% 

0.5 

l.OOE+17 

475 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

1.9671E-02 

8.5910E-01 

1.3959E-02 

7.5000E-01 

1.8612E-02 

82.60% 

13.96% 

0.5 

l.OOE+17 

500 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

1.9671E-02 

8.5874E-01 

1.3915E-02 

7.5000E-01 

1.8553E-02 

82.37% 

13.91% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

**ARC: thickness of ARC is dependent on Lambda-O (optimal wavelength)and ARC refractive index —> arc_depth = Lambda-O/ (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaAs solar cell (N on P on BSF-P+)(AM0) 


Divs*: 

10 

Points: 

4399 

Triangles: 

8528 


Low Doping 





















N-t 

N-d 

P-t 

P-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.5 

1.00E+17 

5 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3510E-02 

8.6683E-01 

1.6939E-02 

7.6000E-01 

2.2883E-02 

83.12% 

12.40% 

0.5 

l.OOE+17 

10 

1.00E+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3429E-02 

8.6441E-01 

1.6842E-02 

7.5999E-01 

2.2160E-02 

83.16% 

12.33% 

0.5 

1.00E+17 

15 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6539E-01 

1.6907E-02 

7.7000E-01 

2.1957E-02 

83.41% 

12.38% 

0.5 

1.00E+17 

20 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6418E-01 

1.6829E-02 

7.6000E-01 

2.2143E-02 

83.14% 

12.32% 

0.5 

l.OOE+17 

25 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6528E-01 

1.6991E-02 

7.6992E-01 

2.2069E-02 

83.84% 

12.44% 

0.5 

1.00E+17 

30 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6577E-01 

1.7037E-02 

7.7000E-01 

2.2125E-02 

84.02% 

12.47% 

0.5 

1.Q0E+17 

35 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6645E-01 

1.7100E-02 

7.7000E-01 

2.2074E-02 

84.25% 

12.52% 

0.5 

l.OOE+17 

40 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3423E-02 

8.6592E-01 

1.7018E-02 

7.7000E-01 

2.2102E-02 

83.91% 

12.46% 

0.5 

1.00E+17 

45 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3424E-02 

8.6635E-01 

1.7063E-02 

7.6999E-01 

2.2159E-02 

84.08% 

12.49% 

0.5 

l.OOE+17 

50 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6543E-01 

1.6935E-02 

7.7000E-01 

2.1994E-02 

83.55% 

12.40% 

0.5 

l.OOE+17 

60 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3404E-02 

8.6483E-01 

1.6842E-02 

7.5999E-01 

2.2161E-02 

83.21% 

12.33% 

0.5 

l.OOE+17 

70 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6527E-01 

1.6902E-02 

7.5999E-01 

2.2240E-02 

83.39% 

12.37% 

0.5 

l.OOE+17 

80 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6455E-01 

1.6815E-02 

7.5999E-01 

2.2126E-02 

83.04% 

12.31% 

0.5 

l.OOE+17 

90 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6490E-01 

1.6856E-02 

7.5999E-01 

2.2179E-02 

83.20% 

12.34% 

0.5 

l.OOE+17 

100 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6426E-01 

1.6783E-02 

7.6000E-01 

2.2083E-02 

82.91% 

12.29% 

0.5 

l.OOE+17 

125 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3406E-02 

8.6407E-01 

1.6755E-02 

7.6000E-01 

2.2047E-02 

82.85% 

12.27% 

0.5 

l.OOE+17 

150 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6412E-01 

1.6786E-02 

7.5999E-01 

2.2087E-02 

82.93% 

12.29% 

0.5 

l.OOE+17 

175 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6399E-01 

1.6743E-02 

7.6000E-01 

2.2030E-02 

82.73% 

12.26% 

0.5 

l.OOE+17 

200 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3410E-02 

8.6390E-01 

1.6741E-02 

7.6000E-01 

2.2028E-02 

82.78% 

12.26% 

0.5 

l.OOE+17 

225 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3426E-02 

8.6390E-01 

1.6728E-02 

7.6000E-01 

2.2010E-02 

82.66% 

12.25% 

0.5 

l.OOE+17 

250 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6391E-01 

1.6750E-02 

7.5999E-01 

2.2040E-02 

82.78% 

12.26% 

0.5 

l.OOE+17 

275 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6376E-01 

1.6710E-02 

7.6000E-01 

2.1986E-02 

82.59% 

12.23% 

0.5 

l.OOE+17 

300 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6411E-01 

1.6787E-02 

7.6000E-01 

2.2188E-02 

82.94% 

12.29% 

0.5 

l.OOE+17 

325 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3423E-02 

8.6378E-01 

1.6722E-02 

7.6000E-01 

2.2003E-02 

82.65% 

12.24% 

0.5 

l.OOE+17 

350 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3425E-02 

8.6370E-01 

1.6695E-02 

7.5999E-01 

2.1967E-02 

82.51% 

12.22% 

0.5 

l.OOE+17 

375 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6416E-01 

1.6789E-02 

7.5999E-01 

2.2091E-02 

82.95% 

12.29% 

0.5 

l.OOE+17 

400 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3412E-02 

8.6382E-01 

1.6721E-02 

7.6000E-01 

2.2001E-02 

82.68% 

12.24% 

0.5 

l.OOE+17 

425 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3424E-02 

8.6366E-01 

1.6690E-02 

7.5997E-01 

2.1961E-02 

82.50% 

12.22% 

0.5 

l.OOE+17 

450 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3423E-02 

8.6446E-01 

1.6830E-02 

7.6000E-01 

2.2145E-02 

83.12% 

12.32% 

0.5 

l.OOE+17 

475 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3422E-02 

8.6412E-01 

1.6774E-02 

7.5999E-01 

2.2072E-02 

82.88% 

12.28% 

0.5 

l.OOE+17 

500 

l.OOE+16 

0.5 

l.OOE+17 

Si02 

0.6000 

2.3421E-02 

8.6385E-01 

1.6722E-02 

7.6000E-01 

2.2003E-02 

82.65% 

12.24% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs’' number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index -> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaAs solar cell (N on P on BSF-P+)(AM1.5) 


Divs*: 

10 

Points: 

4399 

Triangles: 

8528 


High Doping (Davenport, 2004) 


















N-t 

N-d 

P-t 

P-d 

BSF-t 

BSF-d 

ARC*’ 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.2 

1.50E+18 

5 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6476E-02 

9.4821E-01 

1.3383E-02 

8.5000E-01 

1.5745E-02 

85.67% 

13.38% 

0.2 

1.50E+18 

10 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6483E-02 

9.4822E-01 

1.3390E-02 

8.5000E-01 

1.5753E-02 

85.67% 

13.39% 

0.2 

1.50E+18 

15 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6481E-02 

9.4922E-01 

1.3455E-02 

8.5993E-01 

1.5647E-02 

86.01% 

13.46% 

0.2 

1.50E+18 

20 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6481E-02 

9.4821E-01 

1.3387E-02 

8.4999E-01 

1.5749E-02 

85.66% 

13.39% 

0.2 

1.50E+18 

25 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6490E-02 

9.4751E-01 

1.3350E-02 

8.5000E-01 

1.5705E-02 

85.44% 

13.35% 

0.2 

1.50E+18 

30 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6486E-02 

9.4714E-01 

1.3333E-02 

8.4999E-01 

1.5686E-02 

85.39% 

13.33% 

0.2 

1.50E+18 

35 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6497E-02 

9.4707E-01 

1.3341E-02 

8.5000E-01 

1.5695E-02 

85.39% 

13.34% 

0.2 

1.50E+18 

40 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6486E-02 

9.4707E-01 

1.3342E-02 

8.5000E-01 

1.5697E-02 

85.45% 

13.34% 

0.2 

1.50E+18 

45 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6492E-02 

9.4724E-01 

1.3357E-02 

8.5000E-01 

1.5714E-02 

85.50% 

13.36% 

0.2 

1.50E+18 

50 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6483E-02 

9.4730E-01 

1.3359E-02 

8.4999E-01 

1.5717E-02 

85.56% 

13.36% 

0.2 

1.50E+18 

60 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6478E-02 

9.4753E-01 

1.3371E-02 

8.4999E-01 

1.5730E-02 

85.64% 

13.37% 

0.2 

1.50E+18 

70 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6491E-02 

9.4795E-01 

1.3398E-02 

8.5000E-01 

1.5762E-02 

85.70% 

13.40% 

0.2 

1.50E+18 

80 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6480E-02 

9.4803E-01 

1.3395E-02 

8.5000E-01 

1.5759E-02 

85.74% 

13.39% 

0.2 

1.50E+18 

90 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6494E-02 

9.4837E-01 

1.3416E-02 

8.5000E-01 

1.5793E-02 

85.77% 

13.42% 

0.2 

1.50E+18 

100 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6482E-02 

9.4828E-01 

1.3406E-02 

8.5000E-01 

1.5772E-02 

85.77% 

13.41% 

0.2 

1.50E+18 

125 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6482E-02 

9.4845E-01 

1.3410E-02 

8.5000E-01 

1.5777E-02 

85.79% 

13.41% 

0.2 

1.50E+18 

150 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6483E-02 

9.4845E-01 

1.3409E-02 

8.5000E-01 

1.5775E-02 

85.77% 

13.41% 

0.2 

1.50E+18 

175 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6495E-02 

9.4886E-01 

1.3431E-02 

8.5000E-01 

1.5801E-02 

85.82% 

13.43% 

0.2 

1.50E+18 

200 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6485E-02 

9.4865E-01 

1.3416E-02 

8.5000E-01 

1.5784E-02 

85.79% 

13.42% 

0.2 

1.50E+18 

225 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6486E-02 

9.4837E-01 

1.3401E-02 

8.5000E-01 

1.5766E-02 

85.71% 

13.40% 

0.2 

1.50E+18 

250 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6486E-02 

9.4865E-01 

1.3415E-02 

8.5000E-01 

1.5782E-02 

85.78% 

13.41% 

0.2 

1.50E+18 

275 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6495E-02 

9.4891E-01 

1.3431E-02 

8.5000E-01 

1.5802E-02 

85.81% 

13.43% 

0.2 

1.50E+18 

300 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6487E-02 

9.4851E-01 

1.3408E-02 

8.4999E-01 

1.5774E-02 

85.74% 

13.41% 

0.2 

1.50E+18 

325 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6488E-02 

9.4873E-01 

1.3419E-02 

8.5000E-01 

1.5787E-02 

85.78% 

13.42% 

0.2 

1.50E+18 

350 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6497E-02 

9.4894E-01 

1.3433E-02 

8.4999E-01 

1.5803E-02 

85.81% 

13.43% 

0.2 

1.50E+18 

375 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6487E-02 

9.4846E-01 

1.3405E-02 

8.5000E-01 

1.5771E-02 

85.72% 

13.41% 

0.2 

1.50E+18 

400 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6486E-02 

9.4864E-01 

1.3413E-02 

8.5000E-01 

1.5780E-02 

85.76% 

13.41% 

0.2 

1.50E+18 

425 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6493E-02 

9.4882E-01 

1.3425E-02 

8.4999E-01 

1.5794E-02 

85.79% 

13.42% 

0.2 

1.50E+18 

450 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6488E-02 

9.4830E-01 

1.3395E-02 

8.5000E-01 

1.5759E-02 

85.67% 

13.40% 

0.2 

1.50E+18 

475 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6488E-02 

9.4845E-01 

1.3404E-02 

8.4999E-01 

1.5770E-02 

85.72% 

13.40% 

0.2 

1.50E+18 

500 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

1.6487E-02 

9.4860E-01 

1.3411E-02 

8.5000E-01 

1.5778E-02 

85.75% 

13.41% 


‘Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

“ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index —> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaAs solar cell (N on P on BSF-P+)(AM0) 


Divs*: 

10 

Points: 

4399 

Triangles: 

8528 


High Doping (Davenport, 2004) 



















N-t 

N-d 

P-t 

P-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.2 

1.50E+18 

5 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0081E-02 

9.5338E-01 

1.6434E-02 

8.6000E-01 

1.9110E-02 

85.84% 

12.03% 

0.2 

1.50E+18 

10 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0081E-02 

9.5338E-01 

1.6435E-02 

8.6000E-01 

1.9110E-02 

85.84% 

12.03% 

0.2 

1.50E+18 

15 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0072E-02 

9.5413E-01 

1.6508E-02 

8.5999E-01 

1.9196E-02 

86.20% 

12.09% 

0.2 

1.50E+18 

20 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0075E-02 

9.5337E-01 

1.6429E-02 

8.6000E-01 

1.9103E-02 

85.84% 

12.03% 

0.2 

1.50E+18 

25 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0089E-02 

9.5283E-01 

1.6384E-02 

8.5997E-01 

1.9052E-02 

85.60% 

11.99% 

0.2 

1.50E+18 

30 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0084E-02 

9.5253E-01 

1.6362E-02 

8.5995E-01 

1.9027E-02 

85.53% 

11.98% 

0.2 

1.50E+18 

35 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0090E-02 

9.5246E-01 

1.6366E-02 

8.5999E-01 

1.9030E-02 

85.53% 

11.98% 

0.2 

1.50E+18 

40 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0084E-02 

9.5245E-01 

1.6372E-02 

8.5996E-01 

1.9038E-02 

85.59% 

11.99% 

0.2 

1.50E+18 

45 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0061E-02 

9.5254E-01 

1.6365E-02 

8.6000E-01 

1.9029E-02 

85.64% 

11.98% 

0.2 

1.50E+18 

50 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0080E-02 

9.5262E-01 

1.6393E-02 

8.6000E-01 

1.9062E-02 

85.70% 

12.00% 

0.2 

1.50E+18 

60 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0073E-02 

9.5280E-01 

1.6406E-02 

8.6000E-01 

1.9077E-02 

85.78% 

12.01% 

0.2 

1.50E+18 

70 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0090E-02 

9.5312E-01 

1.6441E-02 

8.5999E-01 

1.9117E-02 

85.86% 

12.04% 

0.2 

1.50E+18 

o 

00 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0077E-02 

9.5318E-01 

1.6438E-02 

8.6000E-01 

1.9114E-02 

85.90% 

12.03% 

0.2 

1.50E+18 

90 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0093E-02 

9.5345E-01 

1.6464E-02 

8.6000E-01 

1.9145E-02 

85.94% 

12.05% 

0.2 

1.50E+18 

100 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0080E-02 

9.5339E-01 

1.6453E-02 

8.6000E-01 

1.9131E-02 

85.94% 

12.04% 

0.2 

1.50E+18 

125 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0082E-02 

9.5353E-01 

1.6460E-02 

8.6000E-01 

1.9140E-02 

85.96% 

12.05% 

0.2 

1.50E+18 

150 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0081E-02 

9.5352E-01 

1.6456E-02 

8.6000E-01 

1.9135E-02 

85.94% 

12.05% 

0.2 

1.50E+18 

175 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0094E-02 

9.5384E-01 

1.6485E-02 

8.6000E-01 

1.9168E-02 

86.01% 

12.07% 

0.2 

1.50E+18 

200 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0075E-02 

9.5367E-01 

1.6460E-02 

8.6000E-01 

1.9139E-02 

85.97% 

12.05% 

0.2 

1.50E+18 

225 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0084E-02 

9.5346E-01 

1.6447E-02 

8.6000E-01 

1.9125E-02 

85.89% 

12.04% 

0.2 

1.50E+18 

250 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0086E-02 

9.5368E-01 

1.6466E-02 

8.6000E-01 

1.9147E-02 

85.96% 

12.05% 

0.2 

1.50E+18 

275 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0094E-02 

9.5388E-01 

1.6485E-02 

8.6000E-01 

1.9168E-02 

86.00% 

12.07% 

0.2 

1.50E+18 

300 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0O85E-O2 

9.5357E-01 

1.6456E-02 

8.6000E-01 

1.9134E-02 

85.92% 

12.05% 

0.2 

1.50E+18 

325 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0O89E-O2 

9.5375E-01 

1.6471E-02 

8.5999E-01 

1.9153E-02 

85.97% 

12.06% 

0.2 

1.50E+18 

350 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0096E-02 

9.5391E-01 

1.6486E-02 

8.6000E-01 

1.9170E-02 

86.00% 

12.07% 

0.2 

1.50E+18 

375 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0O86E-O2 

9.5354E-01 

1.6452E-02 

8.6000E-01 

1.9130E-02 

85.90% 

12.04% 

0.2 

1.50E+18 

400 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0078E-02 

9.5367E-01 

1.6457E-02 

8.5999E-01 

1.9136E-02 

85.95% 

12.05% 

0.2 

1.50E+18 

425 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0084E-02 

9.5381E-01 

1.6470E-02 

8.6000E-01 

1.9151E-02 

85.98% 

12.06% 

0.2 

1.50E+18 

450 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0087E-02 

9.5341E-01 

1.6440E-02 

8.6000E-01 

1.9116E-02 

85.84% 

12.04% 

0.2 

1.50E+18 

475 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0086E-02 

9.5353E-01 

1.6451E-02 

8.6000E-01 

1.9129E-02 

85.89% 

12.04% 

0.2 

1.50E+18 

500 

7.00E+17 

0.5 

8.00E+18 

Si02 

0.6000 

2.0086E-02 

9.5365E-01 

1.6460E-02 

8.5998E-01 

1.9140E-02 

85.93% 

12.05% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index -> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaSb solar cell (P on N on BSF-N+)(AM1.5) 


Divs*: 

10 

Points: 

4346 

Triangles: 

8424 


Low Doping 





















P-t 

P-d 

N-t 

N-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.5 

1.00E+17 

5 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

3.6966E-02 

3.2791E-01 

8.8920E-03 

2.7000E-01 

3.2933E-02 

73.36% 

8.89% 

0.5 

1.00E+17 

10 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7013E-02 

3.2810E-01 

8.9103E-03 

2.7000E-01 

3.3001E-02 

73.37% 

8.91% 

0.5 

1.00E+17 

15 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7047E-02 

3.2825E-01 

8.9237E-03 

2.7000E-01 

3.3051E-02 

73.38% 

8.92% 

0.5 

1.00E+17 

20 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7071E-02 

3.2835E-01 

8.9334E-03 

2.7000E-01 

3.3087E-02 

73.39% 

8.93% 

0.5 

l.OOE+17 

25 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7090E-02 

3.2843E-01 

8.9408E-03 

2.7000E-01 

3.3115E-02 

73.40% 

8.94% 

0.5 

1.00E+17 

30 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7103E-02 

3.2849E-01 

8.9461E-03 

2.7000E-01 

3.3134E-02 

73.40% 

8.95% 

0.5 

1.00E+17 

35 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7115E-02 

3.2854E-01 

8.9509E-03 

2.7000E-01 

3.3152E-02 

73.41% 

8.95% 

0.5 

1.00E+17 

40 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7125E-02 

3.2858E-01 

8.9548E-03 

2.7000E-01 

3.3166E-02 

73.41% 

8.95% 

0.5 

1.00E+17 

45 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7131E-02 

3.2862E-01 

8.9576E-03 

2.7000E-01 

3.3176E-02 

73.41% 

8.96% 

0.5 

1.00E+17 

50 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7140E-02 

3.2865E-01 

8.9608E-03 

2.7000E-01 

3.3188E-02 

73.41% 

8.96% 

0.5 

1.00E+17 

60 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7149E-02 

3.2870E-01 

8.9648E-03 

2.7000E-01 

3.3203E-02 

73.42% 

8.96% 

0.5 

1.00E+17 

70 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7161E-02 

3.2874E-01 

8.9691E-03 

2.7000E-01 

3.3219E-02 

73.42% 

8.97% 

0.5 

1.00E+17 

80 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7167E-02 

3.2877E-01 

8.9717E-03 

2.7000E-01 

3.3229E-02 

73.42% 

8.97% 

0.5 

1.00E+17 

90 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7174E-02 

3.2880E-01 

8.9742E-03 

2.7000E-01 

3.3238E-02 

73.42% 

8.97% 

0.5 

1.00E+17 

100 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7175E-02 

3.2881E-01 

8.9750E-03 

2.7000E-01 

3.3241E-02 

73.42% 

8.97% 

0.5 

1.00E+17 

125 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7186E-02 

3.2885E-01 

8.9792E-03 

2.7000E-01 

3.3256E-02 

73.43% 

8.98% 

0.5 

1.00E+17 

150 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7195E-02 

3.2888E-01 

8.9823E-03 

2.7000E-01 

3.3268E-02 

73.43% 

8.98% 

0.5 

1.00E+17 

175 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7201E-02 

3.2890E-01 

8.9842E-03 

2.7000E-01 

3.3275E-02 

73.43% 

8.98% 

0.5 

1.00E+17 

200 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7202E-02 

3.2892E-01 

8.9851E-03 

2.7000E-01 

3.3278E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

225 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7208E-02 

3.2893E-01 

8.9869E-03 

2.7000E-01 

3.3285E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

250 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7201E-02 

3.2899E-01 

8.9855E-03 

2.7000E-01 

3.3280E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

275 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7213E-02 

3.2895E-01 

8.9887E-03 

2.7000E-01 

3.3291E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

300 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7214E-02 

3.2896E-01 

8.9893E-03 

2.7000E-01 

3.3294E-02 

73.43% 

8.99% 

0.5 

1.00E+17 

325 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7215E-02 

3.2897E-01 

8.9896E-03 

2.7000E-01 

3.3295E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

350 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7217E-02 

3.2897E-01 

8.9905E-03 

2.7000E-01 

3.3298E-02 

73.43% 

8.99% 

0.5 

1.0OE+17 

375 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7218E-02 

3.2898E-01 

8.9908E-03 

2.7000E-01 

3.3299E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

400 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7215E-02 

3.2898E-01 

8.9901E-03 

2.7000E-01 

3.3297E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

425 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7218E-02 

3.2898E-01 

8.9910E-03 

2.7000E-01 

3.3300E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

450 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7221E-02 

3.2899E-01 

8.9918E-03 

2.7000E-01 

3.3303E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

475 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7222E-02 

3.2899E-01 

8.9921E-03 

2.7000E-01 

3.3304E-02 

73.43% 

8.99% 

0.5 

l.OOE+17 

500 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

3.7218E-02 

3.2899E-01 

8.9910E-03 

2.7000E-01 

3.3300E-02 

73.43% 

8.99% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength (and ARC refractive index —> arc_depth = Lambda-O/ (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaSb solar cell (P on N on BSF-N+)(AM0) 


Divs*: 

10 

Points: 

4346 

Triangles: 

8424 


Low Doping 




















P-t 

P-d 

N-t 

N-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

imax 

FF 

Eff 

0.5 

1.00E+17 

5 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.6652E-02 

3.3910E-01 

1.4197E-02 

2.8000E-01 

5.0704E-02 

73.90% 

10.39% 

0.5 

1.00E+17 

10 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.6834E-02 

3.3936E-01 

1.4257E-02 

2.8000E-01 

5.0917E-02 

73.92% 

10.44% 

0.5 

1.00E+17 

15 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.6921E-02 

3.3953E-01 

1.4288E-02 

2.8000E-01 

5.1028E-02 

73.93% 

10.46% 

0.5 

1.00E+17 

20 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.6984E-02 

3.3965E-01 

1.4310E-02 

2.8000E-01 

5.1107E-02 

73.94% 

10.48% 

0.5 

1.00E+17 

25 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7039E-02 

3.3975E-01 

1.4329E-02 

2.7999E-01 

5.1177E-02 

73.94% 

10.49% 

0.5 

1.00E+17 

30 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7073E-02 

3.3982E-01 

1.4341E-02 

2.7998E-01 

5.1222E-02 

73.95% 

10.50% 

0.5 

1.00E+17 

35 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7107E-02 

3.3988E-01 

1.4353E-02 

2.7999E-01 

5.1265E-02 

73.95% 

10.51% 

0.5 

1.00E+17 

40 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7130E-02 

3.3993E-01 

1.4361E-02 

2.8000E-01 

5.1291E-02 

73.95% 

10.51% 

0.5 

1.00E+17 

45 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7152E-02 

3.3998E-01 

1.4370E-02 

2.8000E-01 

5.1320E-02 

73.95% 

10.52% 

0.5 

1.00E+17 

50 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7169E-02 

3.4001E-01 

1.4376E-02 

2.8000E-01 

5.1341E-02 

73.96% 

10.52% 

0.5 

1.00E+17 

60 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7197E-02 

3.4005E-01 

1.4386E-02 

2.8000E-01 

5.1378E-02 

73.96% 

10.53% 

0.5 

1.00E+17 

70 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7224E-02 

3.4008E-01 

1.4395E-02 

2.8000E-01 

5.1410E-02 

73.97% 

10.54% 

0.5 

1.00E+17 

80 

3.50E+17 

0.5 

1.00E+18 

Si02 

0.6000 

5.7243E-02 

3.4010E-01 

1.4402E-02 

2.7999E-01 

5.1436E-02 

73.97% 

10.54% 

0.5 

1.00E+17 

90 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7258E-02 

3.4012E-01 

1.4407E-02 

2.7999E-01 

5.1454E-02 

73.98% 

10.55% 

0.5 

1.00E+17 

100 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7269E-02 

3.4014E-01 

1.4411E-02 

2.8000E-01 

5.1467E-02 

73.98% 

10.55% 

0.5 

1.00E+17 

125 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7290E-02 

3.4017E-01 

1.4418E-02 

2.7999E-01 

5.1495E-02 

73.98% 

10.56% 

0.5 

1.00E+17 

150 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7311E-02 

3.4019E-01 

1.4426E-02 

2.8000E-01 

5.1520E-02 

73.99% 

10.56% 

0.5 

1.00E+17 

175 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7324E-02 

3.4021E-01 

1.4430E-02 

2.7999E-01 

5.1537E-02 

73.99% 

10.56% 

0.5 

1.00E+17 

200 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7333E-02 

3.4022E-01 

1.4433E-02 

2.8000E-01 

5.1546E-02 

73.99% 

10.57% 

0.5 

1.00E+17 

225 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7342E-02 

3.4023E-01 

1.4436E-02 

2.8000E-01 

5.1557E-02 

73.99% 

10.57% 

0.5 

1.00E+17 

250 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7341E-02 

3.4024E-01 

1.4436E-02 

2.7999E-01 

5.1559E-02 

73.99% 

10.57% 

0.5 

1.00E+17 

275 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7354E-02 

3.4025E-01 

1.4440E-02 

2.8000E-01 

5.1571E-02 

74.00% 

10.57% 

0.5 

1.00E+17 

300 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7358E-02 

3.4026E-01 

1.4442E-02 

2.8000E-01 

5.1577E-02 

74.00% 

10.57% 

0.5 

1.00E+17 

325 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7362E-02 

3.4026E-01 

1.4443E-02 

2.8000E-01 

5.1582E-02 

74.00% 

10.57% 

0.5 

1.00E+17 

350 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7366E-02 

3.4026E-01 

1.4444E-02 

2.8000E-01 

5.1586E-02 

74.00% 

10.57% 

0.5 

1.00E+17 

375 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7369E-02 

3.4027E-01 

1.4445E-02 

2.8000E-01 

5.1590E-02 

74.00% 

10.57% 

0.5 

1.00E+17 

400 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7371E-02 

3.4027E-01 

1.4446E-02 

2.8000E-01 

5.1592E-02 

74.00% 

10.58% 

0.5 

1.00E+17 

425 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7373E-02 

3.4027E-01 

1.4446E-02 

2.7999E-01 

5.1596E-02 

74.00% 

10.58% 

0.5 

1.00E+17 

450 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7375E-02 

3.4028E-01 

1.4447E-02 

2.8000E-01 

5.1597E-02 

74.00% 

10.58% 

0.5 

1.00E+17 

475 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7377E-02 

3.4028E-01 

1.4448E-02 

2.8000E-01 

5.1600E-02 

74.00% 

10.58% 

0.5 

l.OOE+17 

500 

3.50E+17 

0.5 

l.OOE+18 

Si02 

0.6000 

5.7379E-02 

3.4028E-01 

1.4448E-02 

2.8000E-01 

5.1601E-02 

73.99% 

10.58% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models there 
are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelengthjand ARC refractive index -> arc_depth = Lambda-O/ (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaSb solar cell (P on N on BSF-N+)(AM1.5) 


Divs*: 

10 

Points: 

4346 

Triangles: 

8424 


High Doping (Davenport, 2004) 


















P-t 

P-d 

N-t 

N-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

isc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.37 

9.80E+191 5 11.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8148E-02 

5.0560E-01 

1.5536E-02 

4.3000E-01 

3.6131E-02 

80.55% 

15.54% 

0.37 

9.80E+19 10 1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8217E-02 

5.0607E-01 

1.5579E-02 

4.3000E-01 

3.6231E-02 

80.55% 

15.58% 

0.37 

9.80E+19 15 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8256E-02 

5.0639E-01 

1.5605E-02 

4.3000E-01 

3.6291E-02 

80.55% 

15.61% 

0.37 

9.80E+19 20 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8283E-02 

5.0667E-01 

1.5625E-02 

4.3000E-01 

3.6337E-02 

80.55% 

15.62% 

0.37 

9.8CE+19 25 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8315E-02 

5.0689E-01 

1.5644E-02 

4.3000E-01 

3.6383E-02 

80.55% 

15.64% 

0.37 

9.80E+19 30 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8332E-02 

5.0706E-01 

1.5656E-02 

4.3000E-01 

3.6410E-02 

80.55% 

15.66% 

0.37 

9.80E+19 

35 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8353E-02 

5.0720E-01 

1.5669E-02 

4.3000E-01 

3.6440E-02 

80.55% 

15.67% 

0.37 

9.80E+19 

40 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8366E-02 

5.0732E-01 

1.5678E-02 

4.3000E-01 

3.6460E-02 

80.55% 

15.68% 

0.37 

9.80E+19 45 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8378E-02 

5.0741E-01 

1.5685E-02 

4.3000E-01 

3.6477E-02 

80.55% 

15.69% 

0.37 

9.80E+19 50 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8387E-02 

5.0750E-01 

1.5692E-02 

4.3000E-01 

3.6492E-02 

80.55% 

15.69% 

0.37 

9.8CE+19 60 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8403E-02 

5.0763E-01 

1.5702E-02 

4.3000E-01 

3.6516E-02 

80.54% 

15.70% 

0.37 

9.8CE+19 70 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8415E-02 

5.0774E-01 

1.5710E-02 

4.3000E-01 

3.6534E-02 

80.54% 

15.71% 

0.37 

9.80E+19 80 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8424E-02 

5.0783E-01 

1.5716E-02 

4.3000E-01 

3.6549E-02 

80.54% 

15.72% 

0.37 

9.8CE+19 90 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8425E-02 

5.0789E-01 

1.5718E-02 

4.3000E-01 

3.6554E-02 

80.54% 

15.72% 

0.37 

9.8CE+19 100 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8439E-02 

5.0795E-01 

1.5725E-02 

4.3000E-01 

3.6571E-02 

80.54% 

15.73% 

0.37 

9.80E+19 125 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8445E-02 

5.0806E-01 

1.5731E-02 

4.3000E-01 

3.6584E-02 

80.54% 

15.73% 

0.37 

9.8CE+19 150 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8460E-02 

5.0815E-01 

1.5739E-02 

4.3000E-01 

3.6603E-02 

80.54% 

15.74% 

0.37 

9.8CE+19 175 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8468E-02 

5.0825E-01 

1.5746E-02 

4.3000E-01 

3.6618E-02 

80.53% 

15.75% 

0.37 

9.80E+19 200 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8468E-02 

5.0825E-01 

1.5746E-02 

4.3000E-01 

3.6618E-02 

80.53% 

15.75% 

0.37 

9.80E+19 225 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8476E-02 

5.0829E-01 

1.5750E-02 

4.3000E-01 

3.6627E-02 

80.53% 

15.75% 

0.37 

9.8CE+19 250 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8479E-02 

5.0832E-01 

1.5752E-02 

4.3000E-01 

3.6632E-02 

80.53% 

15.75% 

0.37 

9.80E+19 275 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8479E-02 

5.0834E-01 

1.5752E-02 

4.3000E-01 

3.6634E-02 

80.53% 

15.75% 

0.37 

9.8CE+19 300 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8484E-02 

5.0837E-01 

1.5755E-02 

4.3000E-01 

3.6640E-02 

80.53% 

15.76% 

0.37 

9.8CE+19 325 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8486E-02 

5.0839E-01 

1.5757E-02 

4.3000E-01 

3.6643E-02 

80.53% 

15.76% 

0.37 

9.80E+19 350 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8488E-02 

5.0840E-01 

1.5758E-02 

4.3000E-01 

3.6646E-02 

80.53% 

15.76% 

0.37 

9.8CE+19 375 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8483E-02 

5.0841E-01 

1.5756E-02 

4.3000E-01 

3.6642E-02 

80.53% 

15.76% 

0.37 

9.8CE+19 400 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8491E-02 

5.0843E-01 

1.5760E-02 

4.3000E-01 

3.6650E-02 

80.53% 

15.76% 

0.37 

9.80E+19 425 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8492E-02 

5.0844E-01 

1.5760E-02 

4.3000E-01 

3.6652E-02 

80.53% 

15.76% 

0.37 

9.8CE+19 450 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8493E-02 

5.0845E-01 

1.5761E-02 

4.3000E-01 

3.6653E-02 

80.53% 

15.76% 

0.37 

9.8CE+19 475 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8494E-02 

5.0846E-01 

1.5762E-02 

4.3000E-01 

3.6655E-02 

80.53% 

15.76% 

0.37 

9.80E+191 500 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

3.8489E-02 

5.0846E-01 

1.5760E-02 

4.3000E-01 

3.6651E-02 

80.53% 

15.76% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs’' number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index -> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for GaSb solar cell (P on N on BSF-N+)(AM0) 


Divs*: 

10 

Points: 

4346 

Triangles: 

8424 


High Doping 

Davenport, 2004) 
















P-t 

P-d 

N-t 

N-d 

BSF-t 

BSF-d 

ARC* 

Lambda-0 

Jsc 

Voc 

Pmax 

Vmax 

Imax 

FF 

Eff 

0.37 

9.80E+19 

5 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8293E-02 

5.1661E-01 

2.4343E-02 

4.4000E-01 

5.5324E-02 

80.83% 

17.82% 

0.37 

9.80E+19 

10 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8439E-02 

5.1713E-01 

2.4427E-02 

4.4O00E-01 

5.5516E-02 

80.83% 

17.88% 

0.37 

9.80E+19 

15 

1.00E+20 

15 

1.00E+20 

Si 02 

0.6000 

5.8531E-02 

5.1747E-01 

2.4481E-02 

4.4000E-01 

5 5638E-02 

80.83% 

17.92% 

0.37 

9.80E+19 

20 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8610E-02 

5.1778E-01 

2.4528E-02 

4.4000E-01 

5.5744E-02 

80.82% 

17.96% 

0.37 

9.80E+19 

25 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8657E-02 

5.1801E-01 

2.4557E-02 

4.4000E-01 

5.5812E-02 

80.82% 

17.98% 

0.37 

9.80E+19 

30 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8704E-02 

5.1819E-01 

2.4585E-02 

4.4000E-01 

5.5874E-02 

80.82% 

18.00% 

0.37 

9.80E+19 

35 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8759E-02 

5.1835E-01 

2.4615E-02 

4.4000E-01 

5.5942E-02 

80.82% 

18.02% 

0.37 

9.80E+19 

40 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8790E-02 

5.1848E-01 

2.4633E-02 

4.4000E-01 

5.5984E-02 

80.81% 

18.03% 

0.37 

9.80E+19 

45 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8817E-02 

5.1858E-01 

2.4648E-02 

4.4000E-01 

5.6019E-02 

80.81% 

18.04% 

0.37 

9.80E+19 

50 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8839E-02 

5.1868E-01 

2.4662E-02 

4.4000E-01 

5.6049E-02 

80.81% 

18.05% 

0.37 

9.80E+19 

60 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8876E-02 

5.1882E-01 

2.4683E-02 

4.4000E-01 

5.6098E-02 

80.81% 

18.07% 

0.37 

9.80E+19 

70 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8904E-02 

5.1894E-01 

2.4699E-02 

4.4000E-01 

5.6135E-02 

80.80% 

18.08% 

0.37 

9.80E+19 

80 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8926E-02 

5.1903E-01 

2.4713E-02 

4.4O00E-01 

5.6165E-02 

80.80% 

18.09% 

0.37 

9.80E+19 

90 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8928E-02 

5.1910E-01 

2.4716E-02 

4.4000E-01 

5.6173E-02 

80.80% 

18.09% 

0.37 

9.80E+19 

100 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8960E-02 

5.1917E-01 

2.4733E-02 

4.4000E-01 

5.6210E-02 

80.80% 

18.11% 

0.37 

9.80E+19 

125 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.8989E-02 

5.1930E-01 

2.4750E-02 

4.4000E-01 

5.6249E-02 

80.79% 

18.12% 

0.37 

9.80E+19 

150 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9010E-02 

5.1938E-01 

2.4762E-02 

4.4000E-01 

5.6277E-02 

80.79% 

18.13% 

0.37 

9.80E+19 

175 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9026E-02 

5.1945E-01 

2.4771E-02 

4.4000E-01 

5.6297E-02 

80.79% 

18.13% 

0.37 

9.80E+19 

200 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9018E-02 

5.1949E-01 

2.4769E-02 

4.4000E-01 

5.6294E-02 

80.79% 

18.13% 

0.37 

9.80E+19 

225 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9047E-02 

5.1954E-01 

2.4783E-02 

4.4000E-01 

5.6326E-02 

80.79% 

18.14% 

0.37 

9.80E+19 

250 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9055E-02 

5.1957E-01 

2.4788E-02 

4.4O0OE-O1 

5.6300E-02 

80.79% 

18.15% 

0.37 

9.80E+19 

275 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9061E-02 

5.1960E-01 

2.4792E-02 

4.4000E-01 

5.6344E-02 

80.79% 

18.15% 

0.37 

9.80E+19 

300 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9067E-02 

5.1962E-01 

2.4795E-02 

4.4000E-01 

5.6352E-02 

80.78% 

18.15% 

0.37 

9.80E+19 

325 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9071E-02 

5.1964E-01 

2.4798E-02 

4.4000E-01 

5.6358E-02 

80.78% 

18.15% 

0.37 

9.80E+19 

350 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9075E-02 

5.1966E-01 

2.4800E-02 

4.4000E-01 

5.6363E-02 

80.78% 

18.16% 

0.37 

9.80E+19 

375 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9062E-02 

5.1967E-01 

2.4795E-02 

4.4000E-01 

5.6351E-02 

80.78% 

18.15% 

0.37 

9.80E+19 

400 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9082E-02 

5.1969E-01 

2.4804E-02 

4.4000E-01 

5.6372E-02 

80.78% 

18.16% 

0.37 

9.80E+19 

425 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9085E-02 

5.1970E-01 

2.4805E-02 

4.4000E-01 

5.6376E-02 

80.78% 

18.16% 

0.37 

9.80E+19 

450 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9087E-02 

5.1971E-01 

2.4807E-02 

4.4000E-01 

5.6379E-02 

80.78% 

18.16% 

0.37 

9.80E+19 

475 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9089E-02 

5.1972E-01 

2.4808E-02 

4.4000E-01 

5.6382E-02 

80.78% 

18.16% 

0.37 

9.80E+19 

500 

1.00E+20 

15 

1.00E+20 

Si02 

0.6000 

5.9070E-02 

5.1972E-01 

2.4800E-02 

4.4000E-01 

5.6364E-02 

80.78% 

18.16% 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the single cell models 
there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index -> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for Tandem GaAs (N-P-P+) GaSb (P-N-N+) solar cell (AMO) 
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74.32/ t*4z 

1.14001*00 4,41141-03 7 44051-03 

5.41/ 

74.40 

114*41-02 

1.55/ 

205 55 

0,)7 4.007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

100 

2 41417-02 0.44407-01 1 74337-02 7,70007-01 2 )2047-02 03.11/ 13.13/ 

7.00147-03 4 43357-01 2 57577-03 3.40007-01 4 40447-03 

74.32/ 1.18/ 

1 14001*04 4.40441-03 7.44111-03 

5.41/ 

74,41 

1145*1-02 

t.53/ 

20501 

0,)7 4.007*14 

50 

1.007*20 

15 

1 007*20 

S.02 

too 

2 44057-02 0.44417-01 1 74347-02 7.70007-01 2 32437-02 03.10/ 13.13/ 

4.44137-03 4 43317-01 2 57177-03 3.40007-01 4 54427-03 

74.32/ 1.11/ 

1 14001*00 4.50421-43 7.44021-03 

5.40/ 

74.40 

1 14541-02 

1.53/ 

20507 

0 )7 4 007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 44057-02 0 44447-01 1 74527-02 7 44117-01 2 33207-02 03.17/ 13.14/ 

4 44007-03 4.43217-01 2 54157-03 3.40007-01 4 5*547-03 

74.31/ t**z 

1 15011*00 4.5*501-03 7.43*41-03 

5.50/ 

74.31 

1 14431-02 

t.54/ 

205 20 

0)7 4 007*14 

50 

1.007*20 

15 

1007*20 

Si02 

too 

2 44047-02 0.44747-01 1 74427-02 7 70007-01 2 33017-02 13.12/ 13.13/ 

4.4*147-03 4.43247-01 2 54527-03 3.40007-01 4 57757-03 

74.31/ 1.**/ 

1.14001*00 4.57751-03 7.42001-03 

5.54/ 

74.30 

11*531-02 

1.53/ 

205.07 

0)7 4.007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 41077-02 0.44517-01 1 74107-02 7.40007-01 2 )5447-02 03.05/ 13.11/ 

4.47247-03 4.43227-01 2 541*7-03 3.40007-01 4 54*77-03 

74.31/ 1.»lz 

1 15001*00 4.54171-03 7 55401-03 

5.53/ 

75.54 

1 17521-02 

1.40/ 

204.72 

0 )7 4.007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 44047-02 0,44707-01 1 74717-02 7.70007-01 2 33)47-02 73.25/ 13.14/ 

4 44757-02 4 43217-01 2 55477-03 3.40007-01 4 54347-03 

74.31/ 1*7/ 

1.14001*00 4 54341-03 7.41)51-03 

5.57/ 

74.14 

1 14421-02 

154/ 

205 31 

0 )7 4 007*14 

50 

1.007*20 

15 

1007*20 

S.02 

too 

2 44057-02 0.44407-01 1 74147-02 7 40007-01 2 35717-02 03.03/ 13.12/ 

4 45157-03 4 43157-01 2 55357-03 3.40007-01 4 54747-03 

74.31/ t*7z 

1.8001*00 4 54741-03 7 52051-03 

5.51/ 

75.20 

117401-02 

1.40/ 

20472 

0 )7 4.007*14 

50 

1.007*20 

15 

1007*20 

S.02 

100 

2 44047-02 0.44737-01 1 74247-02 7.40007-01 2 )5147-02 03.04/ 13,12/ 

4.45437-03 4 43177-01 2 55537-03 3.40007-01 4 55217-03 

74.31/ 1.87/ 

1.15001*04 4 55281-03 7 53401-03 

5.52/ 

75.35 

117531-02 

t.40/ 

204 70 

0)7 4,007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 44057-02 0.44447-01 174327-02 7.40007-01 2 35457-02 03.00/ 13.13/ 

4 45157-03 4 43157-01 2 55)47-03 3.40007-01 4 54737-03 

74.31/ 1.87/ 

1 15001*00 4 54731-03 7.52041-03 

551/ 

75.20 

1 17551-02 

1.41/ 

204*4 

0 )7 4 007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 44037-02 0.44417-01 1 74047-02 7 54447-01 2 )5517-02 02.44/ 13.11/ 

4.4)107-03 4.43017-01 2 54557-03 3.40007-01 4 52447-03 

74.31/ I.Hz 

1 15001*00 4 52401-03 7 50501-03 

5.44/ 

75.04 

1 17331-02 

t.5®z 

204 SO 

0)7 4.007*14 

50 

1.007*20 

15 

1007*20 

Si02 

too 

2 44047-02 0.44777-01 1 74237-02 7 40007-01 2 35037-02 03.03/ 13.12/ 

4 44317-03 4.43127-01 2 55027-03 3.40007-01 4.53*47-03 

74.31/ 1.87/ 

1.8001*00 4.53141-03 7 51071-03 

5.50/ 

75.20 

117471-02 

1.40/ 

204.7) 

0)7 4.007*14 

50 

1.007*20 

15 

1007*20 

Si02 

too 

2 44057-02 0.44407-01 1 74247-02 7.40007-01 2 35007-02 03.05/ 13.12/ 

4 43447-03 4.43117-01 2 54177-03 3.40007-01 4 53527-03 

74.31/ 1.87/ 

1.8001*00 4 53521-43 7 5 841-03 

5.50/ 

75.8 

117411-02 

*.40/ 

20475 

0)7 4.007*14 

50 

1.007*20 

15 

1007*20 

Si02 

too 

2 41147-02 0.44507-01 1 71057-02 7 54447-01 2 35)37-02 02.44/ 13.04/ 

4.43347-03 4 43047-01 2 54447-03 3.40007-01 4 52477-03 

74.31/ t**Z 

1.8001*00 4.52071-03 7.50011-03 

5.50/ 

75 00 

117241-02 

1.51/ 

20431 

0.)7 4.007*14 

50 

1.007*20 

15 

1 007*20 

S.02 

too 

2 44047-02 0.44447-01 1 74017-02 7.40007-01 2 35547-02 02.44/ 13.10/ 

4.43217-03 4.43017-01 2 54547-03 3.40007-01 4 52*07-03 

74.31/ t*iz 

1.8001*00 4 52101-03 7 50721-03 

5.50/ 

75.07 

117321-02 

1.50/ 

20447 

0.)7 4.007*14 

50 

1.007*20 

15 

1 007*20 

Si02 

too 

2 44047-02 0.44707-01 1 74417-02 7,40007-01 2 34047-02 03.11/ 13.13/ 

4.42417-03 4.43077-01 2 54507-03 3.40007-01 4 52547-03 

74.31/ 1.80/ 

1 8001*00 4 52541-03 7 50451-03 

5.40/ 

75.04 

117511-02 

t.40/ 

204 14 

0)7 4.007*14 

50 

1.007*20 

15 

1007*20 

Si02 

too 

2 44057-02 0.44477-01 1 74117-02 7.40007-01 2 35407-02 02.40/ 13.11/ 

4.42707-03 4.43047-01 2 54)47-03 3.40007-01 4 52247-03 

74.31/ t*tz 

1 15001*00 4 52201-03 7 50131-03 

5.44/ 

75.01 

117351-02 

1.50/ 

204 55 

0 )7 4 007*14 

50 

1007*20 

« 

1 007*20 

Si02 

too 

2 41157-02 0.44547-01 1 71757-02 7 40007-01 2 35207-02 02.04/ 13.04/ 

4 41257-03 4 43017-01 2 53*37-03 3.40007-01 4 50137-03 

74.31/ 1*0/ 

1 8001*00 4 50131-03 7 4*451-03 

541/ 

74*5 

117111-02 

1.57/ 

204 14 


Combined Results 


•Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the 
model (in the single cell models there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index --> arc.depth = lambda-0 / (4*arc_refractivejndex) 
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SILVACO ATLAS Modelling Stats for Tandem GaAs (N-P-P+) GaSb (P-N-N+) solar cell (AMO) 


PamVr 4344 1 .hhWi 1424 


7-U M»l 


Upper Cell (GaAs) 

N-d P-i P-d BSF-i BSF-d ARC’ 



Lower Cell (GaSb) 

P-« P-d N-t N-d BSF-t BSF-d ARC" L-0 


j 1MC2* 

11 **€•:• 
imci* 


m 


I I MCI* 
1MC2* 
IMCI* 
n 14K.2* 
I MCI* 
I MCI* 
I MCI* 
I MCI* 
I MCI* 
| I MCI* 
I I MCI* 
1 I MCI* 


1MCI0 

1*K"24 
I-**€•!* 
1 MCI* 


l MCI* KOI 


MCI* tiOI i 


I MCI* 
1MC24 
IMCI* 
I MCI* 
I MCI* 
1MCI* 
I MCI* 
I MC24 
I MCI* 
1MC2* 
IMCI* 
I MCI* 
1MC2* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
1 MCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
IMCI* 
I MCI* 


m iiEB m 3 B M i 


14K*24 SMI 
l*K>2* SM I 


■**€•10 SMI 


. **€•!« SMI 


I. **€•!* SMI 
IMCI* SMI 
I *♦€•!* SMI 
I.MC»I* SMI 
I-**€•!« SMI 
1MC»I* SMI 
1MC.I* SMI 
IMC»I* SMI 


1*K«24 SMI 
1 **€•!* SMI 


I.MC'I* SMI 


Upper Cell 

Voc Pmax Vmax Imax FF 


no) 31-01 4 14301-04 

I.43721-41 5 0*10104 
I.M4C-41 141)71**) 
4.437K-41 l4342t-*2 


3 5MOC-OI 5.52*41-44 

7.40MC-41 7.4INC-44 
7 oMOE-41 OI47K-43 
7M«C-»I 134*11-02 


7S.74X 

M.3IX 

74.43X 


unit «i 
2 3)231-02 

2 04421-02 

2.45711-02 
2.45721-42 
2 0*421-02 
2*5411*2 
2.43*21-42 
2 *372t-02 
2.43! It'*2 
2*2311-42 


0 255K-41 141*11*2 
I I43U-01 141301**2 
».34«4t-«1 1**031*2 

0*1471-01 144*01*2 

0*2031-01 14*101*2 

♦ 53*21-41 UHX*! 
0.547K-41 147451-42 
o **1*1-01 U47K-*2 

♦ 55551-01 147121-42 

♦*41*101 147301 *2 

*.54431-41 147121-42 
*.54341-*l 144*11*2 
0*»OK-OI HHIt*I 

♦ 53*11-41 15MU-42 


-tHt 


..IV 


♦***1*1 I774K-42 
7J44K-4I 1H4K-42 
I.***K-41 1425K-42 
1 40071-01 1*3441-42 
UMOt-OI 1*5451-42 
I5HHJ1 10,101-42 
HWMI 104721-02 
I.IMKtl t.*3*4C-*2 
MMOE-OI 1.44321-42 
M*MC-«1 M44K-42 
I.4MK-4I 1,4431c**! 
I4MM-41 M1371-42 
14*401-01 1.17*41-41 
1.444K-41 1745K-02 


45.34X 

14.44/ 

IS.**/ 


I.»IMC'*2 

24H11-02 

IMR4I 
2.44*71-42 
2 4M31-02 

1. »I13C-*I 

2. »**•€■*! 
201101-02 
I.«t*IC-»2 


2.»ll4C-*2 

2.41431-42 

2.01171-02 

2.41421-02 

i*w4t-»: 

I.*MK'*2 

2.41111-42 

1. *mc*2 

2. »**l€-*2 
2411*1-02 

201111-02 
2 00*01-02 
2 41131-02 
2 41241-02 
2411*1-02 
24MT1-02 
2.41171-42 
2.4II1C-*! 
2 01*01-02 
2 41121-42 


♦ 55*21-41 U 7)11-42 

0*)*0C-01 1*4441*2 

0 533K-01 114411-*! 

♦ 53421-41 114551-*! 

4.53411**1 1*4*11-42 
4 S34IC-01 14452€*«2 
♦.534C-41 1*4551-42 

4.S417C**! 1**301-42 
4,534*C**1 U4*2t-42 
4.I2MC-41 1*40*1*2 
*525*1-41 M3131-42 
4.5244C**! t*,*0t-*2 
0*2411-01 M3*21-42 
4.52SIC**! 1*3071-42 
4,I245C**1 1*4141-42 
4*2131*1 1*4211-42 
0*31*1-01 1*4*11-42 
4.5322C**! 1*4*01-42 
4*341101 1*41*1*2 
♦.*,421-41 L447IC-*! 
*53551-41 144741 *2 
0*3**101 144771 *2 
*.5317t-41 1*5051-42 

♦ *3711-01 144171*2 

4.S354C-01 1*4*11-42 
4*3711-01 1*4171-42 

4,53411-01 14*0*1 *2 
4.S344C-41 1*47*1-42 
4.537IC-41 I444!t-*I 
4.5344C-41 1**071-42 
4.5357C-*! 1*4721-42 
4S37U-01 144IIC-*! 

0*31*1-41 1444IC-*! 

♦ 5,441-41 14440C-*2 
4.5354C-01 1.447H-*! 
4 *3411-01 144111 *2 


I.4M0C-01 1 4374C-02 H.HX 
M44K-41 1.412 3t-02 *S.T4z 
l.4M*C-*1 14134C-0I 15.14/ 
I.5444C-*! 1.412 4C-02 1*14/ 
I4M0E-01 1013K-42 15.14/ 

M44K-41 1.41341-02 15.15/ 
4.5444C-01 1.422K-02 MUX 
I4M0E-0I 1413*02 *5.14/ 
l.4«*0C-*1 1.4*75C-*I *5.4*/ 
«4M*C-*I 1.044K-42 *5.53/ 
*.*44*1-01 1.405IC02 *5.53/ 
I.4M0C-M 100**02 15.5*Z 

I.4M0C-01 14*441 *2 *5.44/ 

•.4M0C-01 1001*1-02 *5.7*/ 

I.SOOIC-OI M1431-42 *5.7*/ 
l,*44tt-*l 1.4142C-OI ts.nx 
I *4441-01 1.4140C-02 *5.44/ 

4.4M0C-01 14144C-02 15.44/ 

I.4MK-M 1.4153C-02 15.44/ 

1.54441- 41 1.4142C-02 15.44/ 

I.544IC-41 1 41*41-02 15.44/ 

*.4M«e-*1 1.4142C-4I 14.01/ 

I.4M4C-41 1.417H-02 15.47/ 

I.5444C-41 1.4144C-02 I5.«4/ 
4.5444C-41 1.417K-0I 15.44/ 
•■4M4C-41 141421-02 M.M/ 
1.5444C-41 1.4154C-02 *5.42/ 
I.4M4C-41 1.41741-02 *5.47/ 
**4441-41 14144C-0I M.M/ 
*.4*MC-*1 1.4154C-02 tS.Mx 
I.4M0C-01 1.414 5C-02 15.45/ 
1.5441C-41 1.41141-0! *5.41/ 

«.4M0C-»1 14144C-02 *5.15/ 

1.54441- 01 1.41*31-02 *5.»4/ 
*4M*t-»1 14U5C-OI *5.43/ 


LowerCell 

Voc Pmax Vmax Imax FF EfM 


4.4234C-41 *5*47t-*3 

4.4174C-01 «.3234t-*3 
4.4440C-01 I.M7H-4I 
4.1445C-01 7.554«€-*3 


4.2M0C-0I 2.03 731-02 

42*MC-*1 1.41 MC-42 
42*MC«1 14O42C-02 
4.2M4C-41 1.7441C-42 


:a: 

I74,3t-.I 4 1744C-0, .4.141-0, 4M.0C-*. 1.7*2.t-OI 7,,*/ 


144331-02 
14M2C-02 
142371-02 
U2MC-42 
141131-02 
1.41* It-*2 
14135C-42 
1*1121-02 
1.4*7«-*2 
11MK-02 


4171H-41 4.41*4C-*3 
4.1*17C-01 4-5744€-*3 
4.1*021-01 * 34141-*, 
4.14741-01 4 24401-«3 
4.I47C-01 4 2*27t-*3 
4.1447t-01 42714t-*3 
4.14431 *1 42*171-43 
4 I454C-41 4 2*20t-*3 
4.14541-01 4 243It-*3 
4 14441-01 4 2240t-*3 
4.143 IC-01 *14*lt-*3 
4.I424C-4I 4,1727t-*3 
4I4IK-01 412*7t-*3 
4 *3431-41 4.*t07t-«3 


4 IMOt-OI 14I7H-42 
4.WMC-M 1.40471-02 
44M4C-M 1.*5*4t-*2 
4.*M4C-*I 1 *3*11-02 
4. IMOt-OI 1.5324C-42 
4M*4t-4t 15247t-42 

4.1M*t-*l 1.52721-42 
41MM 41 1*2441-42 

4-1M4t-*1 15227t-42 
4*M0t-01 151IH-42 

4 IMOt-OI 1 51I4C-42 

4 4M4C-41 154551-42 
4.*M4t-41 1 444U-42 

4-1M4t-*1 14I3K-0I 


74.44/ 

74.44/ 

74.44/ 

74.44/ 


7.1547C-43 
7.135*1-43 
7.11*31-42 
7.H34C-43 
7.4*551-43 
7.4*4 31-43 
7.47411-43 
7.44451-43 
7.MM1-43 
7.45HC-43 


7.43451-43 
7.MM1-43 
7.43141-03 
7.*24*t-*3 
7.41751-43 
7.41**1-43 


4.4**41-43 

4.447*1-43 

4.44751-03 

4.45KC-43 

4.45431-43 

4.45151-43 

4.443*t-*3 

4.44311-43 

4.43441-43 

4.431«-*3 

4.43*51-43 

4.424IC-43 

4.427*1-43 


4 44321-41 2.4*441-«3 
4.43471-41 2-44451-*3 
4.431*1-01 2*3111-43 

4.43741- 41 2 42531-«3 

4.43741- 41 241141-0, 
4.43751-41 2.42441-43 
4.434H-41 2 41351 *3 
4.43431-41 2.4*4*l-*3 
4,43511-01 2*0131-03 
4.43541-01 2.54451-*3 
443*41-01 254451-43 
4*3*01-01 2.54251-«3 
4.43441-41 2.54141-*3 
4.43451-41 2.5*741-«3 
4.43441-41 2*tl3l-03 
4 43431-01 2**441 *3 
4.434K-41 2 5*2*1-«3 
4 431*1-41 2.57421-*3 
4.41371-41 2*71*1-03 
4.43331-41 2*7341-03 
4,43241-41 2.54431-*! 
4.432*1-41 2.54*51-«3 
4.41241-41 2*4521-0 3 
4 432K-01 2 5*411*1 
4.432K-41 2**471-03 
4,43171-01 2.55421-«3 
4*3171-01 2.5551I-*! 
4 43151-01 2 5*341-43 
4.43121-*1 2 55021 *3 
4.41121-41 2 *5021-0, 
4.431K-41 2 54*71 *3 
4.43*11-41 2 54551-*1 
4*3011-01 2*4531 *3 
4.43*71-41 2-54541-«3 
4.43*41-41 2 54341-0, 
4,43*51-51 2*4241-03 


3.4M41-41 4.112*1-03 
3.0M4t-*l 4 71*11-01 
1.4***t-*l 4.744*1-43 
3.4*44t-*l 4.73*41-43 
3.4M4C41 4.71*11-41 
3.4M4I-41 4.71441-43 
3.4M4C41 4.74H1-43 
3 4M0I-41 4.4*411-43 
3.4*4*t-*1 4.44441-43 
3.4**4€-*1 4 44*31 *3 
3 4M01-01 4 4*771-43 
3.4*MI-*1 4.44741-43 
3.4M41-41 4.44511-43 
3.4M41-*! 4 434)1-41 
3.4**4€-*1 * *3*71-43 
3.4M4I-*! 4.42711-43 
3.4M4t-*l 4.422*1-43 
3.0M4I-41 4.41341-43 
3.0M41-41 4.41*41-*3 
5.OM0I-41 4.5*041-43 
3.MMI-4I 4.51 IK-43 
3.0M01-41 4*1*01-01 
3.0M01-01 4*7751-03 
3.0M4I-4I 4.54171-41 
3.0M4C-*! * *4341-43 
3.0M41-41 4 55441-43 
3.MM1-0I 45*211-43 
3.0M41-4I 4.54731-43 
3.4MK-M 4 53*01-43 
3.0M4I-4I 4.53IK-43 
3.MM1-41 4*3*21-03 
3.4M4I-41 * *2*01-43 
3.0M41-*! 4 *2441-43 
3.MM1-41 4.52541-43 
3.MMI-41 4.52201-43 
3.0M41-41 5.32*21-43 


74.12/ 
7*43/ 
70.12/ 
70.32/ 
70.32/ 
70.12/ 
70.12/ 
70.32/ 
70.32/ 
70.32/ 
70.12/ 
70.32/ 
70.32/ 
70.12/ 
70.32/ 
70.32z 
70.31/ 
70.11/ 


t*7x 

1.17/ 

1.47/ 

MIX 

MTz 

1.47/ 

1*7/ 

1.17/ 

1.14/ 

Mix 

Mi/ 

1.14/ 

I**/ 


Combined Results 


Dual Junction 


P-i EJf-i W/mT2 


> 5 *2041-04 4.44031-*4 

> 7 011*1-04 0.25*41-*4 

I 0.147*1-03 1H10C-02 

> 1.34*11-02 1 4 0141-02 


12H01-M 

M4Mt*M 
122MCM 
124MC-M 
124MC»M 
12*071.04 
127MCM 
127MC-M 
tITMCM 
127MCM 
L27MCM 
127MI-00 

1 2TMCM 

127M1-M 
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125MCM 
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125MCM 
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125MCM 
12*0*1.00 
125*K«M 
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125MCM 
125M1.M 
125MI.M 
12*MCM 
125M1.M 
125M1.M 
125MCM 
125«0C0O 
1.254K.M 
125Mt« M 
1.25M1.M 
125MCOO 
125M1.M 
1.25MCM 
125MC«M 


4.1I2K-03 14*2*1-*) 
4.7**K-*3 4.474K-43 
4 74*51-03 *.433«-«3 
4,73*41-47 1.41451-43 
4.75143 I 30311-03 
4.700K-43 1 30031-03 
4.7*44t-*3 1 37471 *3 
4.4«4K-*3 I 3*511-03 
4.44OK-43 I.337K-43 
4.44531-43 *.3314t-«3 
4.45771-43 *.322«-*3 
4.447K-43 t 34011-03 
4.44511-43 *.3*721-*3 
4.43431-43 I.202K-43 
4.43471-43 *ao5K-«3 
4.427)1-03 UI3K-43 
4.42211-43 1^7131-03 
4.413K-43 1 24471-03 
4.41*41-43 I2405103 
4 *0001-03 I 24041-0, 
4.5IIK-43 1 23401-03 
4.5I5K-43 I 23221-03 
5.57751-43 122111-03 
4.54371-43 12*4*1-0) 
4*4)41-03 120421-0, 
4 55441-03 I.103K-43 
4.S52C-43 I 10021-0, 
4.54731-43 I.114K-43 
4 53IK-0, 1.17,41-0) 

4.53IK-43 1 173K-0, 
4*3*21-03 t.HOK-43 
4*2*01-03 I.15I41-*) 
4.52441-03 I.157K-43 
4*2*01-03 *.1S7K-*3 
4.522K-43 1.15)41-03 
4.52*21-43 *.15*21-43 


4.43/ 

i.«X 

4.42/ 

4- *1Z 
4*1/ 
4.**/ 

4. MZ 

5- MX 

5. M/ 
5.**/ 


M.42 

*4.7* 

44.33 


*3.77 

•3.55 

*3.37 

*3.32 

13.22 

43.M 

43.47 

*2.43 

42.44 

42.14 

*2.7* 

*2.47 

*2.45 

*2.5* 

•2.35 

42.32 

*222 

*2.45 

*2.44 


Tandem 

P-t Elf-t W/m‘2 


I.7IIO1-0, 4.43/ 

*.45411-43 4.34/ 

121331-42 UK 
14*141-42 11.72/ 


2.44)31-42 

2 44101-02 14.45X 

1141*1-42 13.43/ 

1*4241-42 13.MX 
143421-42 14.44/ 

143*11-42 14.0/ 

144221-42 14.22/ 

143441-42 14.24/ 

1.03441-42 14.11/ 

1433K-*2 14.44/ 
142151-42 14.12/ 
141451-42 44.45/ 

1.01201*02 14.MZ 

1*4*41-42 13.4*/ 

14*351-42 13.74/ 


iMS 


11-42 7 


10*501-02 

1M371-42 

14*2*1-02 

1M331-42 

10*401-02 

1M251-42 

10*201-02 

1M341-02 

14MK-42 

1M201-02 

14*271-42 

14MK-42 

1.MOI-42 


211.11 
2*7.1* 
210 55 
227.53 
221.42 
220.4* 
23*.«7 
224.31 
224.55 
224.45 


1.42711-42 7.52/ 

1*1*21-41 7,47*/ 

14MK-42 7.34/ 

14MK-42 7)1/ 

1.M741-42 7.3*/ 

1.4*111-42 7.34/ 

1M741-42 
1.4*431-42 
14*421-42 
1.M34I-42 
14*241-42 
14*251-42 
14*241-42 7.34/ 

14*231-42 
1M321-42 
1.M35I-42 
14*411-02 
1.4*441-42 
1.4*541-42 
14*441-42 
14*431-42 


1M.57 

1M21 

10*27 


Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in the 


single cell models there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

••ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelength)and ARC refractive index -> arc_depth = Lambda-0 / (4*arc_refractive_index) 
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SILVACO ATLAS Modelling Stats for Tandem GaAs (N-P-P+) GaSb (P-N-N+) solar cell (2000 K Blackbody) 
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‘Divisions: used as a measure of model complexity, a setting for the number of vertical and horizontal divisions for each region in the model (in 
the single cell models there are about 5 X-axis, and 10 Y-axis divisions, each subdivided by the "divs" number). 

“ARC: thickness of ARC is dependent on Lambda-0 (optimal wavelengthjand ARC refractive index --> arc depth=Lambda-0 / (4*arc refractive jndex) 
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APPENDIX B. SINGLE CELL DECKBUILD SCRIPT 


********************************************************* 

# gaas cell test plots using Atlas # 

# Thomas, M: Aug 2015 # 

######################################################### 

go atlas 

############################################# 

## Parameters # 

############################################# 

# Cell Width (urn) 
set cell_width=10 

# Dimensional multipliers (for realistic current and power numbers) 

# Scaling up from set size to 1 square centimeter, 
set full_width_mult=10000/$cell_width 

set length_mult=10000 

set J_mult=$full_width_mult*$length_mult 

# Standard # of Divisions [for easy model scaling] 

# An increase of one unit represents a unit squared increase in model complexity (at 
least) 

# Minimum successfully tested: 3 [only took a couple minutes — 1 and 2 completed, but 
produced strange model results]. 

# Maximum successfully tested: 30 [took about 34 minutes — above that ran out of memory 
or had too many triangles]. 

# [System Config: 8 core 2.4 GHz processor, 12 GB RAM] 
set divs=10 

# Contact (%) 

set contact_percent=8 

# Contact Specifications (cathode = upper contact, anode = lower contact) 
set c_mat="Aluminum" 

set a_mat="Aluminum" 
set c_depth=0.3 
set a_depth=0.3 

set c_width=$cell_width*(0.01*$contact_percent) 

#set c_width=$cell_width*0.99 

set c_x_min=($cell_width/2)-($c_width/2) 

set c_x_max=($cell_width/2)+($c_width/2) 

set a_width=$cell_width*(0.01*$contact_percent) 

#set a_width=$cell_width*0.99 

set a_x_min=($cell_width/2)-($a_width/2) 

set a_x_max=($cell_width/2)+($a_width/2) 

#Anti-Reflective Coating (ARC) thickness 

#Per pveducation.org: optimal AR coating thickness (d) = (wavelength) / (4*refractive 
index) 

# d = Lambda_0 / (4*n) 

# to minimize reflection at that wavelength (0.6 micron in this case). 

#Per refractiveindex.info: n(Si02) = 1.4585, n(Ti02) = 2.6142, n(Zr02) = 2.1248, n(A1203) 
= 1.7557 

# [for testing purposes, n for Air is 1.0002742 

#Per IOPScience "Antireflection coatings for GaAs solar cell applications, Zr02 and A1203 
are 

# components of a highly effective ARC for GaAs solar cells 
set Lambda_0 = 0.6 
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set arc_mat="Si02" 

set arc_ref_index=l.7557 

set arc_depth=$Lambda_0/(4*$arc_ref_index) 

#ARC around contacts 

set c_arc_width=($cell_width-$c_width)/2 
set a_arc_width=($cell_width-$a_width)/2 

#Material for gaps above ARC and next to Contacts 
set gap_mat="Vacuum" 

#GaAs (n, p, and bsf): 

# N-Type layer 
set n_mat="GaAs" 
set n_depth=0.2 

set n_dope=l.5000e+18 

# P-Type layer 
set p_mat="GaAs" 
set p_depth=100 

set p_dope=7.0000e+17 

# BSF Layer 

set bsf_mat="GaAs" 
set bsf_depth=0.5 
set bsf_dope=8.000e+18 

#Upper and Lower cell limits (n, p, and bsf regions) 
set j_depth=0 

set l_top=$j_depth-$n_depth 

set l_bot=$j_depth+$p_depth+$bsf_depth 

############################################# 

#Mesh specification # 

############################################# 


# The "scale_slide" values are used to dynamically adjust the density of the model 

# in the substrate layer, p (for upper cell, n for lower cell) in this case. 

# For very small values of layer depth (e.g. below 1 micron), these values can cause 

# problems with the model, commenting them out here and where they're used in the 

# y.mesh specifications can allow the model to run with fairly similar results. 

set n_p_scale_slide=($n_depth/($p_depth/2)) 
set p_bsf_scale_slide=($bsf_depth/($p_depth/2)) 

# The width mulitplier for the mesh uses a combination of the multipliers for both 

# width and length scaling; while this is not physically accurate, it results in 

# correct values for amperage output (in A/cm A 2). 

mesh width=$J_mult 


#— X divisions 

x.mesh loc=0.0 

x.mesh loc=$c_arc_width 

x.mesh loc=$c_arc_width+$c_width 

x. mesh loc=$cell_width 
#— Y divisions 

y. mesh loc=$l_top-$c_depth 
y.mesh loc=$l_top-$arc_depth 
y.mesh loc=$l_top 

y.mesh loc=$j_depth 

y.mesh loc=$j_depth+$n_p_scale_slide 
y.mesh loc=$j_depth+($p_depth/2) 
y.mesh loc=$l_bot-$p_bsf_scale_slide 
y.mesh loc=$l_bot 
y.mesh loc=$l_bot+$arc_depth 
y.mesh loc=$l_bot+$a_depth 


s=$c_arc_width/$divs 
s=$c_width/$divs 
s=$c_width/$divs 
s=$c_arc_width/$divs 

s=$ c_dept h/$ a r c_dept h 
s=$arc_depth 
s=$n_depth/$divs 
s=$n_depth/$divs 
s=$n_depth/$divs 
s=$p_depth/$divs 
s=$bsf_depth/$divs 
s=$bsf_depth/$divs 
s=$ a_dept h/$ a r c_dept h 
s=$a_depth/$arc_depth 
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##Regions 

region num=l material=$gap_mat x.min=0.0 x.max=$c_x_min y.min=$l_top- 

$c_depth y.max=$l_top-$arc_depth 

region num=2 material=$arc_mat x.min=0.0 x.max=$c_x_min y.min=$l_top- 

$arc_depth y.max=$l_top 

region num=3 material=$c_mat x.min=$c_x_min x.max=$c_x_max y.min=$l_top-$c_depth 
y.max=$l_top 

region num=4 material=$arc_mat x.min=$c_x_max x.max=$cell_width y.min=$l_top- 
$arc_depth y.max=$l_top 

region num=5 material=$gap_mat x.min=$c_x_max x.max=$cell_width y.min=$l_top- 

$c_depth y.max=$l_top-$arc_depth 

region num=6 material=$n_mat x.min=0.0 x.max=$cell_width y.min=$l_top 

y.max=$ j_depth 

region num=7 material=$p_mat x.min=0.0 x.max=$cell_width y.min=$j_depth 

y.max=$ j_depth+$p_depth 

region num=8 material=$bsf_mat x.min=0.0 x.max=$cell_width y.min=$l_bot- 

$bsf_depth y.max=$l_bot 

region num=9 material=$gap_mat x.min=0.0 x.max=$a_x_min y.min=$l_bot+$arc_depth 

y.max=$l_bot+$a_depth 

region num=10 material=$arc_mat x.min=0.0 x.max=$a_x_min y.min=$l_bot 

y.max=$l_bot+$arc_depth 

region num=ll material=$a_mat x.min=$a_x_min x.max=$a_x_max y.min=$l_bot 
y.max=$l_bot+$a_depth 

region num=12 material=$arc_mat x.min=$a_x_max x.max=$cell_width y.min=$l_bot 
y.max=$l_bot+$arc_depth 

region num=13 material=$gap_mat x.min=$a_x_max x.max=$cell_width 
y.min=$l_bot+$arc_depth y.max=$l_bot+$a_depth 

##Electrodes 

electrode num=l name=anode region=3 material=$c_mat x.min=$c_x_min x.max=$c_x_max 
y.min=$l_top-$c_depth y.max=$l_top 

electrode num=2 name=cathode region=ll material=$a_mat x.min=$a_x_min x.max=$a_x_max 
y.min=$l_bot y.max=$l_bot+$a_depth 

##Doping 
#— N 

doping gaussian region=6 n.type conc=$n_dope 
#— P 

doping gaussian region=7 p.type conc=$p_dope 
#— BSF 

doping gaussian region=8 p.type conc=$bsf_dope 

# set contact material to be opaque 
material material=$c_mat imag.index=1000 
material material=$a_mat imag.index=1000 

# :T: Struct out after initial completion [for testing] 

#structure outfile=gaas_test.str 

############################################# 

#Testing # 

############################################# 

# Models 

models conmob fldmob consrh print 

# Method 

method gummel newton trap atrap=0.8 maxtraps=20 climit=le-4 

# set light beam using inbuilt solar spectrum (AMI.5 or AMO) 

#beam num=l x.origin=$cell_width/2 y.origin=$l_top-10 angle=90 AMI.5 back.ref1 
beam num=l x.origin=$cell_width/2 y.origin=$l_top-10 angle=90 AMO back.ref1 

# - 

# test for IV curve: 

# - 
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solve init 
solve previous 

log outfile=gaas_test_iv.log 

solve bl=l beam=l vcathode=-0.01 vstep=0.01 vfinal=1.0 name=cathode 
#tonyplot gaas_test_iv.log -set iv_curve.set 

# :T: Struct out after jsc, voc, iv curve tests to get photogeneration results [for 
testing] 

structure outfile=gaas_test.str 

# - 

# test for spectral response: 

# - 

solve init 
solve previous 

log outfile=gaas_test_sr.log 

solve bl=l beam=l lambda=0.30 wstep=0.01 wfinal=1.2 
#tonyplot gaas_test_sr.log -set spec_resp.set 

#- 

# calculate internal and external quantum efficiency: 

# - 

extract init inf="gaas_test_sr.log" 

extract name="EQE" curve(elect."optical wavelength", -(i."cathode")/elect."source photo 
current") outf="gaas_test_EQE.dat" 

extract name="IQE" curve(elect."optical wavelength", -(i."cathode")/elect."available 
photo current") outf="gaas_test_IQE.dat" 

#tonyplot gaas_test_EQE.dat -overlay gaas_test_IQE.dat -set quant_eff.set 

# - 

# calculate overall cell efficiency and general electrical parameters): 

# - 

# The below is adapted from the SILVACO example solarexl2 

# input variables are appropriately adjusted for the init infile (e.g. cathode vice anode 
—> whatever 

# some variable names are adjusted for clarity as well (e.g. "Pm" —> "Pmax") 

#Power in is based on integrated value of desired spectrum converted from W/m / '2 —> 

W/cm A 2 

#For reference: AMO = 1366 W/m A 2, AMI.5g (global) = 1000 W/m A 2, AMI.5d (direct, for 
concentrators) = 900 W/m A 2 
#set P_In = (900 / 10000) 

#set P_In = (1000 / 10000) 
set P_In = (1366 / 10000) 

extract init infile="gaas_test_iv.log" 

extract name="Jsc" y.val from curve(v."cathode", i."anode") where x.val=0.0 

extract name="Voc" x.val from curve(v."cathode", i."anode") where y.val=0.0 

extract name="P" curve(v."cathode", (v."cathode" * i."anode")) outf="gaas_test_P.dat" 

extract name="Pmax" max(curve(v."cathode", (v."cathode" * i."anode"))) 

extract name="Vmax" x.val from curve(v."cathode", (v."cathode"*i."anode")) where 

y.val=$"Pmax" 

extract name="Imax" $"Pmax"/$"Vmax" 

extract name="FF" $"Pmax"/($"Jsc"*$"Voc")*100 

extract name="Eff" ($Pmax/$P_In)*100 

#Real world GaAs (thin film, 0.9927 cm A 2) stats for comparison (Jsc in mA/cm A 2 —> 

A/cm A 2, Voc in V): 

#Per "Solar cell efficiency tables (Version 45)", Green, et. al; Prog. Photovolt: Res. 

Appl. 2015; 23: 1-9 

set Jsc_realCell = (29.68/1000) 

set Voc_realCell = 1.122 
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set FF_realCell = 0.865*100 
set Eff_realCell = 0.288*100 


extract name 
extract name 
extract name 
extract name 


= " J s c. 

_check" 

$ Jsc 

/ 

$Jsc_realCell 

= "Voc 

_check" 

$Voc 

/ 

$Voc_realCell 

= "FF_i 

check" 

$FF 

/ 

$FF_realCell 

= "Eff 

_check" 

$Ef f 

/ 

$Eff_realCell 


# :T: Struct out after all tests run [for testing] 
♦structure outfile=gaas_test.str 

#tonyplot gaas_test.str 
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APPENDIX C. TANDEM CELL DECKBUILD SCRIPT 


********************************************************* 

# tandem gaas-gasb cell test plots using Atlas # 

# Thomas, M: Aug 2015 # 

######################################################### 

go atlas 

############################################# 

## Parameters # 

############################################# 

# Cell Width (urn) 
set cell_width=10 

# Dimensional multipliers (for realistic current and power numbers) 

# Scaling up from set size to 1 square centimeter, 
set full_width_mult=10000/$cell_width 

set length_mult=10000 

set J_mult=$full_width_mult*$length_mult 

# Standard # of Divisions [for easy model scaling] 

# An increase of one unit represents a unit squared increase in model complexity (at 
least) 

# Minimum successfully tested: 3 [only took a couple minutes — 1 and 2 completed, but 
produced strange model results]. 

# Maximum successfully tested: 30 [took about 34 minutes — above that ran out of memory 
or had too many triangles]. 

# [System Config: 8 core 2.4 GHz processor, 12 GB RAM] 
set divs=10 

# Contact (%) 

set contact_percent=8 

# Contact Specifications (cathode = upper contact, anode = lower contact; u=upper cell, 
l=lower cell) 

set uc_mat="Aluminum" 
set ua_mat="Aluminum" 
set uc_depth=0.5 
set ua_depth=0.5 

set lc_mat="Aluminum" 
set la_mat="Aluminum" 
set lc_depth=0.5 
set la_depth=0.5 

set uc_width=$cell_width*(0.01*$contact_percent) 

#set uc_width=$cell_width*0.99 

set uc_x_min=($cell_width/2)-($uc_width/2) 

set uc_x_max=($cell_width/2)+($uc_width/2) 

set ua_width=$cell_width*(0.01*$contact_percent) 

#set ua_width=$cell_width*0.99 

set ua_x_min=($cell_width/2)-($ua_width/2) 

set ua_x_max=($cell_width/2)+($ua_width/2) 

set lc_width=$cell_width*(0.01*$contact_percent) 

#set lc_width=$cell_width*0.99 

set lc_x_min=($cell_width/2)-($lc_width/2) 

set lc_x_max=($cell_width/2)+($lc_width/2) 

#set la_width=$cell_width*(0.01*$contact_percent) 


83 



set la_width=$cell_width*0.99 

set la_x_min=($cell_width/2)-($la_width/2) 

set la_x_max=($cell_width/2)+($la_width/2) 


(d) = (wavelength) / (4*refractive 


#Anti-Reflective Coating (ARC) thickness 
#Per pveducation.org: optimal AR coating thickness 
index) 

# 

# to minimize reflection at that wavelength (0.6 
#Per refractiveindex.info: n for Si02 is 1.4585, n for Ti02 is 2.6142 

# [for testing purposes, n for Air is 1.0002742 

#After testing Si02, Ti02, and Air, Si02 is an effective ARC for Silicon cells. 


d = 
and 1. 


Lambda_0 / (4*n) 

0 micron in this case). 


#Upper cell ARC 

set u_Lambda_0 = 0.6 

set uarc_mat="Si02" 

set uarc_ref_index=l.4585 

set uarc_depth=$u_Lambda_0/(4*$uarc_ref_index) 


# ARC around contacts 

set uc_arc_width=($cell_width-$uc_width)/2 
set ua_arc_width=($cell_width-$ua_width)/2 


#Lower cell ARC 

#set l_Lambda_0 = (1.0 - $u_Lambda_0) 
set l_Lambda_0 = 1.0 
set larc_mat="Si02" 
set larc_ref_index=l.4585 

set larc_depth=($l_Lambda_0/(4*$larc_ref_index)) 
# ARC around contacts 

set lc_arc_width=($cell_width-$lc_width)/2 
set la_arc_width=($cell_width-$la_width)/2 


#Material for gaps above ARC and next to Contacts 
set gap_mat="Vacuum" 

#Material and gap distance for between cells 
set tgap_mat="Ti02" 
set tgap_depth=0.5 

#P, N, and BSF layers: 

#Upper cell 

#- 

# N-Type layer 
set un_mat="GaAs" 
set un_depth=0.2 

set un_dope=l.5000e+18 

# P-Type layer 
set up_mat="GaAs" 
set up_depth=l.0 

set up_dope=7.0000e+17 

# BSF Layer (P+) 
set ubsf_mat="GaAs" 
set ubsf_depth=0.5 

set ubsf_dope=8.000e+18 

#Lower cell 

# - 

# P-Type layer 
set lp_mat="GaSb" 
set lp_depth=0.37 

set lp_dope=9.8000e+19 

# N-Type layer 
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set ln_mat =,, GaSb n 

set ln_depth=45 

set ln_dope=l.0000e+20 

# BSF Layer (N+) 
set lbsf_mat="GaSb" 
set lbsf_depth=0.37 
set lbsf_dope=l.0000e+20 

#Upper and Lower cell limits (p, n, and bsf regions — for some simpler / shorter math in 
mesh statements) 
set uj_depth=0 

set u_top=$uj_depth-$un_depth 

set u_bot=$uj_depth+$up_depth+$ubsf_depth 

set 1j_depth=$u_bot+$ua_depth+$tgap_depth+$lc_depth+$lp_depth 

set l_top=$lj_depth-$lp_depth 

set l_bot=$lj_depth+$ln_depth+$lbsf_depth 

############################################# 

#Mesh specification # 

############################################# 


# The "scale_slide" values are used to dynamically adjust the density of the model 

# in the substrate layer, p (for upper cell, n for lower cell) in this case. 

# For very small values of layer depth (e.g. below 1 micron), these values can cause 

# problems with the model, commenting them out here and where they're used in the 

# y.mesh specifications can allow the model to run with fairly similar results. 

set un_p_scale_slide=($un_depth/($up_depth/2)) 
set up_bsf_scale_slide=($ubsf_depth/($up_depth/2)) 

set lp_n_scale_slide=($lp_depth/($ln_depth/2)) 
set ln_bsf_scale_slide=($lbsf_depth/($ln_depth/2)) 

# The width mulitplier for the mesh uses a combination of the multipliers for both 

# width and length scaling; while this is not physically accurate, it results in 

# correct values for amperage output (in A/cm A 2). 

mesh width=$J_mult 


#— X divisions 

x.mesh loc=0.0 

x.mesh loc=$uc_arc_width 

x.mesh loc=$uc_arc_width+$uc_width 

x.mesh loc=$cell_width 


s=$uc_arc_width/$divs 
s=$uc_width/$divs 
s=$uc_width/$divs 
s=$uc_arc_width/$divs 


#— Y divisions 

#-Upper Cell 

y.mesh loc=$u_top-$uc_depth 
y.mesh loc=$u_top-$uarc_depth 
y.mesh loc=$u_top 
y.mesh loc=$uj_depth 

y.mesh loc=$uj_depth+$un_p_scale_slide 
y.mesh loc=$uj_depth+($up_depth/2) 
y.mesh loc=$u_bot-$up_bsf_scale_slide 
y.mesh loc=$u_bot 
y.mesh loc=$u_bot+$uarc_depth 
y.mesh loc=$u_bot+$ua_depth 


s=$uc_depth/$uarc_depth 
s=$uarc_depth 
s=$un_depth/$divs 
s=$un_depth/$divs 
s=$un_depth/$divs 
s=$up_depth/$divs 
s=$ubsf_depth/$divs 
s=$ubsf_depth/$divs 
s=$ua_depth/$uarc_depth 
s=$ua_depth/$uarc_depth 


#-Tandem Cell Gap: average between bottom of top and top of bottom cells 

y.mesh loc=(($u_bot+$ua_depth)+($l_top-$lc_depth))/2 s=$tgap_depth/$divs 


#-Lower Cell 

y.mesh loc=$l_top-$lc_depth 
y.mesh loc=$l_top-$larc_depth 
y.mesh loc=$l_top 
y.mesh loc=$lj_depth 


s=$lc_depth/$larc_depth 
s=$larc_depth 
s=$lp_depth/$divs 
s=$lp_depth/$divs 
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y.mesh loc=$lj_depth+$lp_n_scale_slide 
y.mesh loc=$lj_depth+($ln_depth/2) 
y.mesh loc=$l_bot-$ln_bsf_scale_slide 
y.mesh loc=$l_bot 
y.mesh loc=$l_bot+$larc_depth 
y.mesh loc=$l_bot+$la_depth 


s=$lp_depth/$divs 

s=$ln_depth/$divs 

s=$lbsf_depth/$divs 

s=$lbsf_depth/$divs 

s=$la_depth/$larc_depth 

s=$la_depth/$larc_depth 


##Regions 
#—Upper Cell 

region num=l material=$gap_mat x.min=0.0 

$uc_depth y.max=$u_top-$uarc_depth 
region num=2 material=$uarc_mat x.min=0.0 
$uarc_depth y.max=$u_top 

region num=3 material=$uc_mat x.min=$uc_x_min 

$uc_depth y.max=$u_top 

region num=4 material=$uarc_mat x.min=$uc_x_max 
$uarc_depth y.max=$u_top 

region num=5 material=$gap_mat x.min=$uc_x_max 

$uc_depth y.max=$u_top-$uarc_depth 
region num=6 material=$un_mat x.min=0.0 

y.max=$uj_depth 

region num=7 material=$up_mat x.min=0.0 

y.max=$uj_depth+$up_depth 

region num=8 material=$ubsf_mat x.min=0.0 
$ubsf_depth y.max=$u_bot 

region num=9 material=$tgap_mat x.min=0.0 
y.min=$u_bot+$uarc_depth y.max=$u_bot+$ua_depth 
region num=10 material=$uarc_mat x.min=0.0 
y.max=$u_bot+$uarc_depth 

region num=ll material=$ua_mat x.min=$ua_x_min 

y.max=$u_bot+$ua_depth 

region num=12 material=$uarc_mat x.min=$ua_x_max 
y.max=$u_bot+$uarc_depth 

region num=13 material=$tgap_mat x.min=$ua_x_max 
y.min=$u_bot+$uarc_depth y.max=$u_bot+$ua_depth 


x.max=$uc_x_min 
x.max=$uc_x_min 
x.max=$uc_x_max 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$ua_x_min 
x.max=$ua_x_min 
x.max=$ua_x_max 
x.max=$cell_width 
x.max=$cell_width 


#—Tandem Cell Gap 

region num=14 material=$tgap_mat x.min=0.0 x.max=$cell_width 

y.min=$u_bot+$ua_depth y.max=$l_top-$lc_depth 


#—Lower Cell 

region num=15 material=$tgap_mat x.min=0.0 
$lc_depth y.max=$l_top-$larc_depth 
region num=16 material=$larc_mat x.min=0.0 
$larc_depth y.max=$l_top 

region num=17 material=$lc_mat x.min=$lc_x_min 
$lc_depth y.max=$l_top 

region num=18 material=$larc_mat x.min=$lc_x_max 
$larc_depth y.max=$l_top 

region num=19 material=$tgap_mat x.min=$lc_x_max 
$lc_depth y.max=$l_top-$larc_depth 
region num=20 material=$lp_mat x.min=0.0 
y.max=$lj_depth 

region num=21 material=$ln_mat x.min=0.0 
y.max=$lj_depth+$ln_depth 

region num=22 material=$lbsf_mat x.min=0.0 
$lbsf_depth y.max=$l_bot 

region num=23 material=$gap_mat x.min=0.0 
y.min=$l_bot+$larc_depth y.max=$l_bot+$la_depth 
region num=24 material=$larc_mat x.min=0.0 
y.max=$l_bot+$larc_depth 

region num=25 material=$la_mat x.min=$la_x_min 
y.max=$l_bot+$la_depth 

region num=26 material=$larc_mat x.min=$la_x_max 
y.max=$l_bot+$larc_depth 

region num=27 material=$gap_mat x.min=$la_x_max 
y.min=$l_bot+$larc_depth y.max=$l_bot+$la_depth 


x.max=$lc_x_min 
x.max=$lc_x_min 
x.max=$1c_x_max 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$cell_width 
x.max=$la_x_min 
x.max=$la_x_min 
x.max=$1a_x_max 
x.max=$cell_width 
x.max=$cell_width 


y.min=$u_top- 
y.min=$u_top- 
y.min=$u_top- 
y.min=$u_top- 
y.min=$u_top- 
y.min=$u_top 
y.min=$uj_depth 
y.min=$u_bot- 

y.min=$u_bot 
y.min=$u_bot 
y.min=$u_bot 


y.min=$l_top- 
y.min=$l_top- 
y.min=$l_top- 
y.min=$l_top- 
y.min=$l_top- 
y.min=$l_top 
y.min=$lj_depth 
y.min=$l_bot- 

y.min=$l_bot 
y.min=$l_bot 
y.min=$l_bot 
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##Electrodes (actual contact names set here — adjust here if necessary, not in mesh 
specification) 

electrode num=l name=anode region=3 material=$uc_mat x.min=$uc_x_min x.max=$uc_x_max 

y.min=$u_top-$uc_depth y.max=$u_top 

electrode num=2 name=cathode region=ll material=$ua_mat x.min=$ua_x_min x.max=$ua_x_max 
y.min=$u_bot y.max=$u_bot+$ua_depth 

electrode num=3 name=collector region=17 material=$lc_mat x.min=$lc_x_min x.max=$lc_x_max 
y.min=$l_top-$lc_depth y.max=$l_top 

electrode num=4 name=emitter region=25 material=$la_mat x.min=$la_x_min x.max=$la_x_max 
y.min=$l_bot y.max=$l_bot+$la_depth 

##Doping 

#— Upper Cell: GaAs (N on P on P+) 

# -N 

doping gaussian region=6 n.type conc=$un_dope 

# -P 

doping gaussian region=7 p.type conc=$up_dope 

# -BSF 

doping gaussian region=8 p.type conc=$ubsf_dope 

#— Lower Cell: GaSb (P on N on N+) 

# -P 

doping gaussian region=20 p.type conc=$lp_dope 

# -N 

doping gaussian region=21 n.type conc=$ln_dope 

# -BSF 

doping gaussian region=22 n.type conc=$lbsf_dope 

# set contact material to be opaque 
material material=$uc_mat imag.index=1000 
material material=$ua_mat imag.index=1000 
material material=$lc_mat imag.index=1000 
material material=$la_mat imag.index=1000 

# :T: Struct out after initial completion [for testing] 
structure outfile=gaas_gasb_tandem.str 

############################################# 

#Testing # 

############################################# 

# Models 

models conmob fldmob consrh print 

# Method 

method gummel newton trap atrap=0.8 maxtraps=20 climit=le-4 

# set light beam using solar spectrum from external file 

#beam num=l x.origin=$cell_width/2 y.origin=$u_top-100 angle=90 AMI.5 back.ref1 
beam num=l x.origin=$cell_width/2 y.origin=$u_top-100 angle=90 AMO back.ref1 

# - 

# test for IV curve: 

# - 

#—Upper cell 
solve init 
solve previous 

log outfile=gaas_gasb_tandem_iv_u.log 

solve bl=l beam=l vcathode=-0.01 vstep=0.01 vfinal=1.0 name=cathode 

#—Lower cell 
solve init 
solve previous 

log outfile=gaas_gasb_tandem_iv_l.log 
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solve bl=l beam=l vcollector=-0.01 vstep=0.01 vfinal=1.0 name=collector 

#tonyplot gaas_gasb_tandem_iv_u.log -overlay gaas_gasb_tandem_iv_l -set iv_curve.set 

# :T: Struct out after jsc, voc, iv curve tests to get photogeneration results [for 
testing] 

structure outfile=gaas_gasb_tandem.str 

# - 

# test for spectral response: 

# - 

solve init 
solve previous 

log outfile=gaas_gasb_tandem_sr.log 

solve bl=l beam=l lambda=0.30 wstep=0.01 wfinal=2.0 

#tonyplot gaas_gasb_tandem_sr.log -set spec_resp.set 

#- 

# calculate internal and external quantum efficiency: 

# - 

extract init inf="gaas_gasb_tandem_sr.log" 

extract name="EQE_u" curve(elect."optical wavelength", (i."anode")/elect."source photo 
current") outf="gaas_gasb_tandem_EQE_u.dat" 

extract name="IQE_u" curve(elect."optical wavelength", (i."anode")/elect."available photo 
current") outf="gaas_gasb_tandem_IQE_u.dat" 

extract name="EQE_l" curve(elect."optical wavelength", (i."emitter")/elect."source photo 
current") outf="gaas_gasb_tandem_EQE_l.dat" 

extract name="IQE_l" curve(elect."optical wavelength", (i."emitter")/elect."available 
photo current") outf="gaas_gasb_tandem_IQE_l.dat" 

extract name="EQE_t" curve(elect."optical wavelength", (i."anode" + 
i."emitter")/elect."source photo current") outf="gaas_gasb_tandem_EQE_t.dat" 
extract name="IQE_t" curve(elect."optical wavelength", (i."anode" + 
i."emitter")/elect."available photo current") outf="gaas_gasb_tandem_IQE_t.dat" 

#tonyplot gaas_gasb_tandem_EQE_u.dat -overlay gaas_gasb_tandem_IQE_u.dat \ 

# -overlay gaas_gasb_tandem_IQE_l.dat -overlay gaas_gasb_tandem_IQE_l.dat \ 

# -overlay gaas_gasb_tandem_IQE_t.dat -overlay gaas_gasb_tandem_IQE_t.dat -set 
quant_eff.set 

# - 

# calculate overall cell efficiency and general electrical parameters): 

# - 

# The below is adapted from the SILVACO example solarexl2 

# input variables are appropriately adjusted for the init infile (e.g. cathode vice anode 
—> whatever 

# some variable names are adjusted for clarity as well (e.g. "Pm" —> "Pmax") 

#Power in is based on integrated value of desired spectrum converted from W/m / '2 —> 

W/cm / "2 

#For reference: AMO = 1366 W/m A 2, AMI. 5g (global) = 1000 W/m /> '2, AMI. 5d (direct, for 
concentrators) = 900 W/m^2 
#set P_In = (900 / 10000) 

#set P_In = (1000 / 10000) 
set P_In = (1366 / 10000) 


# Calculations for upper cell # 

# - # 

extract init infile="gaas_gasb_tandem_iv_u.log" 

extract name="uJsc" y.val from curve(v."cathode", i."anode") where x.val=0.0 
extract name="uVoc" x.val from curve(v."cathode", i."anode") where y.val=0.0 
extract name="uP" curve(v."cathode", (v."cathode" * i."anode")) 

outf="gaas_gasb_tandem_uP.dat" 

extract name="uPmax" max(curve(v."cathode", (v."cathode" * i."anode"))) 
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extract name="uVmax" x.val from curve(v."cathode", (v."cathode"*i."anode")) where 
y.val=$"uPmax" 

extract name="ulmax" $"uPmax"/$"uVmax" 

extract name="uFF" $"uPmax"/($"uJsc"*$"uVoc")*100 

extract name =,, uEff" ($uPmax/$P_In) *100 

#Real world GaAs (thin film, 0.9927 cm A 2) stats for comparison (Jsc in mA/cm /v 2 —> 
A/cm /> '2, Voc in V) : 

#Per "Solar cell efficiency tables (Version 45)", Green, et. al; Prog. Photovolt: Res. 

Appl. 2015; 23: 1-9 

set Jsc_realCellu = (29.68/1000) 

set Voc_realCellu = 1.122 

set FF_realCellu = 0.865*100 

set Eff_realCellu = 0.288*100 

extract name ="Jsc_check" $uJsc / $Jsc_realCellu 
extract name ="Voc_check" $uVoc / $Voc_realCellu 
extract name ="FF_check" $uFF / $FF_realCellu 
extract name ="Eff_check" $uEff / $Eff_realCellu 


#-# 

# Calculations for lower cell # 

# - # 


extract init infile="gaas_gasb_tandem_iv_l.log" 

extract name="lJsc" y.val from curve(v."collector", i."emitter") where x.val=0.0 
extract name="lVoc" x.val from curve(v."collector", i."emitter") where y.val=0.0 
extract name="lP" curve(v."emitter", (v."collector" * i."emitter")) 
outf="gaas_gasb_tandem_lP.dat" 

extract name="lPmax" max(curve(v."collector", (v."collector" * i."emitter"))) 

extract name="lVmax" x.val from curve(v."collector", (v."collector"*i."emitter")) where 

y.val=$"lPmax" 

extract name="Umax" $"lPmax"/$"lVmax" 

extract name="lFF" $"lPmax"/($"1Jsc"*$"lVoc")*100 

extract name="lEff" ($lPmax/$P_In)*100 

#Real world GaSb stats for comparison (Jsc in A/cm A 2, Voc in V): 

#Per "Fundamental Characterization Studies of GaSb Solar Cells," Fraas et. al; Boeing 

High Technology Center, Seattle, WA. 

set Jsc_realCelll = (5.5/1000) 

set Voc_realCelll = 0.491 

set FF_realCelll = 0.813*100 

set Eff_realCelll = 0.07*100 

extract name ="Jsc_check" $lJsc / $Jsc_realCelll 
extract name ="Voc_check" $lVoc / $Voc_realCelll 
extract name ="FF_check" $1FF / $FF_realCelll 
extract name ="Eff_check" $lEff / $Eff_realCelll 


#-# 

# Results for combined cell performance, as dual-junction or # 

# tandem design, were calculated in an external application # 

# (e.g. Microsoft Excel) based on upper and lower cell outputs. # 

# - # 


# :T: Struct out after all tests run [for testing] 
#structure outfile=gaas_gasb_tandem.str 

#tonyplot gaas_gasb_tandem.str 

quit 
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